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INTRODUCTION 

Improved  molecular  imaging  tools  to  further  research  and  development  and  to  translate  into  the  clinic 
are  needed  to  improve  breast  cancer  disease  diagnosis  and  treatment.  Radioluminescence  Imaging  (RLI) 
/  Tomography  (RLT)  have  been  developed  in  this  grant  to  accomplish  these  goals.  RLI  utilizes  optical 
light  emitted  from  either  radionuclides  or  radio-luminescent  biologically-compatible  phosphor 
nanoparticles  (RLNPs)  to  produce  images  of  molecular  contrasts;  these  molecular  contrasts  are  intended 
to  provide  expansive  disease-specific  information.  This  has  particularly  relevance  for  breast  cancer 
because  of  the  limited  specificity  of  current  breast  cancer  imaging  tools.  Because  radiation  sources 
commonly  used  in  medicine,  such  as  radionuclides,  x-ray  sources,  and  linear  accelerators  provide  the 
excitation,  RLI  has  the  potential  to  function  online  during  conventional  imaging  procedures  or  surgical 
procedures. 

The  work  completed  in  this  fellowship  completed  the  following  aims: 

Aim  1:  Select  and  characterize  the  optimal  nano-scintillators  for  X-ray  imaging  in  tissue 

Aim  2:  Develop  an  adaptable  optical  imaging  system,  functional  under  X-ray  irradiation,  to 

enable  the  detection  of  emitted  light  from  nano-scintillators.  Develop  the  reconstruction  imaging 

software  to  perform  the  X-ray/optical  reconstruction. 

Aim  3:  Analyze  the  performance  characteristics  of  Nano-scintillator  imaging  with  tissue- 

simulating  phantoms.  Compare  to  Optical  Fluorescence  Imaging. 

Aim  4:  Determine  the  feasibility  of  injectable  nano-scintillators  in  vivo. 

This  grant  has  not  only  provided  for  extensive  investigation  into  specific  techniques,  parameters, 
and  limits  of  RLI  as  a  tool  to  accomplish  these  goals,  but  has  also  led  to  the  grantee’s  intensive  study  in 
molecular  targeting  of  cancer,  medical  physics  and  radiation  therapy,  and  an  opportunity  to  mentor 
undergraduates,  graduate  students,  and  postdoctoral  fellows  in  RLI  and  molecular  imaging.  Aims  1-4 
were  completed.  These  efforts  resulted  in  13  journal  publications,  including  5  first  author  (2  under 
consideration),  12  conference  abstracts,  and  8  oral  presentations. 

BODY 

1.1.  Nanophosphor  characterization 

To  meet  the  milestones  of  Aiml,  the  feasibility  of  X-ray  luminescence  with  Gadolinium  oxy sulfide 
particles  was  investigated.  We  examined  the  practical  aspects  of  this  new  modality,  including  phosphor 
concentration,  light  emission  linearity,  detector  damage,  and  spectral  emission  characteristics.  Finally, 
the  contrast  produced  by  these  phosphors  was  compared  to  that  of  X-ray  fluoroscopy. 

Figure  1  shows  the  characteristics  of  gadolinium  (GOS)  and  lanthanum  (LOS)  oxysulfide 
phosphors  doped  with  terbium  (green),  europium  (red),  thulium  (NIR),  and  Praseodymiun  (NIR).  These 
phosphors  were  studied,  due  to  their  anticipated  brightness.  Figure  la  demonstrates  the  GOS:Tb 
emission.  The  light  emission  was  imaged  in  a  clinical  X-ray  scanner  with  a  cooled  CCD  camera  and  a 
spectrophotometer;  dose  measurements  were  determined  with  a  calibrated  dosimeter;  the  setup  is  shown 
in  Figure  lb. 

Dose  to  the  CCD  camera  in  the  chosen  imaging  geometry  was  measured  at  less  than 
0.02cGy/sec.  Emitted  light  was  found  to  be  linear  with  dose  (R  =  1)  and  concentration  (R  =1),  as 
shown  in  Figure  Ic.  Phosphor  emission  peaks  were  all  less  than  20nm  full-width  at  half-maximum,  and 
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less  than  3nm  full- width  at  half-maximum  for  Eu  and  Tb,  as  was  expected  from  lanthanide  dopants,  as 
showin  in  Figure  Id.  The  minimum  practical  concentration  necessary  to  detect  luminescent  phosphors 
was  dependent  on  dose;  it  was  estimated  that  sub-picomolar  concentrations  are  detectable  at  the  surface 
of  the  tissue  with  typical  mammographic  doses,  with  the  minimum  detectable  concentration  increasing 
with  depth  and  decreasing  with  dose.  In  a  reflection  geometry.  X-ray  luminescence  had  nearly  a  430- 
fold  greater  contrast  to  background  than  X-ray  contrast. 

The  outcome  of  this  aim  was  a  successful  demonstration  and  feasibility  assessment  of  this 
modality  using  Gadolinium  Oxysulfide.  However,  due  to  difficulties  in  making  these  phosphors 
biocompatible,  the  phosphor  platform  was  switched  to  BaYF4,  which  may  be  more  easily  bioconjugated 
and  tagged  to  molecular  markers.  Indeed,  in  a  preliminary  study,  these  phosphor  were  targeted  to  the 
Folate  receptor  (commonly  expressed  in  breast  cancer),  and  uptaken  by  live  cells  (data  not  shown). 


Wavelength  (nm) 


Figure  1  (a)  Image  of  a  vial  of  RFNP  mixed  with  agarose  excited  by  x-ray.  (b)  Excitation  setup  for 
image  in  (a),  (c)  Concentration  vs.  photon  emission  linearity,  (d)  Normalized  emission  spectra  from 
FaOS:  doped  with  Terbium  (Tb),  Praseodymium  (Pr),  Thulium  (Tm),  and  Europium  (Eu). 


1.2:  Nanophosphor  characterization 

In  performing  these  tasks,  the  investigator  was  been  exposed  to  the  field  of  molecular  imaging,  a  new 
direction  for  this  PI.  This  research  education  was  aided  with  participation  in  BioE222:  Molecular 
Imaging,  which  brought  together  the  molecular  imaging  faculty  at  Stanford  to  teach  aspects  in  the 
hardware,  chemistry,  and  biology  of  molecular  imaging.  In  addition,  the  PI  was  exposed  to  nanoparticle 
fabrication,  including  the  processes  in  making  nanoparticles  stable  in  human  serum  with  low  toxicity. 
Also,  the  PI  gained  knowledge  in  molecular  targets,  and  the  advantages  and  disadvantages  of  targeting 
to  peptides,  hormones,  antibidies,  affabodies,  and  other  targeting  agents.  This  program  has  been  aided 
by  working  at  benchside  with  these  materials  scientists,  biologists,  and  nuclear  imaging  experts. 

2.1:  Imaging  Hardware  Development 
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Aim  2  extended  the  utility  of  Radioluminescent  Imaging  with  the  construction  of  a  fully  functioning 
small-animal  imaging  scanner,  shown  in  Figure  2,  which  has  full  spectral  and  imaging  capabilities  to 
image  small-animals  injected  with  nanophosphors  with  high  resolution.  Figure  2  shows  the  small- 
animal  imaging  scanner,  with  a  rod- lens  installed.  The  system  was  created  by  repurposing  a  retired 
commercial  optical  imaging  scanner,  by  adding  a  more  sensitive  camera  (EM-CCD,  Pro-EM,  Princeton 
Instruments),  a  controllable  ambient  light  source,  and  filters  appropriate  for  RENP  phosphors. 


In  addition,  we  improved  the  sensitivity  of 
Radioluminescent  Imaging  by  incorporating  leaded 
glass  in  front  of  the  camera,  a  chiller,  and  a  filter 
wheel.  The  3/8”  thick  glass  sheet  was  custom 
machined  (Ray-bar,  Inc)  to  attach  to  our  REI  imager, 
and  attenuates  the  X-rays  by  >99.999%  before  they  hit 
the  CCD  and  cause  increased  noise  in  the  pixels.  Due 
to  this  improvement,  we  have  seen  a  dramatic  decrease 
in  the  stochastic  X-ray  noise  in  the  CCD;  this  has 
enabled  integration  times  of  over  20  seconds.  The 
liquid  chiller  for  the  CCD  camera  brings  the  CCD 
temperature  down  to  -85C. 

We  demonstrated  the  improvements  to  the  REI 
system  by  imaging  an  organic  phosphor  (PDOT)  and 
water  (H20)  with  the  old  setup  and  the  new  setup.  We 
placed  the  samples  in  our  imaging  box  and  irradiated 
them  with  50  and  80  kVp  x-rays  using  a  superficial  X- 
ray  irradiator.  Eor  the  old  samples,  we  shut  off  the 
chiller  and  removed  the  X-ray  glass.  Using  ideal  gain 
and  integration  settings,  signal  for  the  PDOTs  and 
H20  were  nearly  indistinguishable  from  the  noise  of  about  200  counts,  due  to  the  high  X-ray  noise  from 
the  incident  irradiation,  and  the  dark  counts  from  thermal  noise  on  the  CCD.  After  connecting  the 
chiller  and  installing  the  leaded  glass  in  front  of  the  CCD  lens,  we  then  repeated  the  experiment  with 
ideal  settings.  As  shown  in  Eigure  3,  signal  for  the  phosphor  was  increased  by  a  factor  of  8  over  the  old 
setup. 


Figure  2:  Small-animal  imaging  box 
fabricated  to  provide  proof-of-concept  and 
aid  in  technique  and  nanoparticle 
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Figure  3:  RLT  system  with  improved  sensitivity,  (a)  3/8”  thick  leaded  glass  prevents  X-ray 
noise  from  hitting  the  CCD,  and  enables  longer  integration  times.  The  filter  wheel  enables 
spectral  selectivity,  (b)  Intensity  emitted  from  phosphors  (PDOT)  and  water.  The  increase 
in  intensity  of  the  (new)  vs  (old)  samples  demonstrates  the  longer  integration  times  possible. 


2.2:  Imaging  Software  Reconstruction  Development 

In  addition,  this  aim  resulted  in  the  development  of  a  novel  algorithm  which  incorporated  a  diffuse 
optical  photon  propagation  model  into  the  reconstruction  algorithm  to  recover  unresolved  dimensions  in 
an  X-ray  limited  angle  (LA)  geometry.  This  sophisticated  technique  enables  such  applications  as  image- 
guided  surgery,  where  the  ability  to  resolve  lesions  at  depths  of  several  centimetres,  which  can  be  the 
key  to  successful  resection. 

The  goal  of  X-ray  Luminescence  Tomography  is  to  determine  the  phosphor  distribution.  This 
spatial  distribution  can  be  determined  by  minimizing  the  difference  between  the  measured  photon  flux 
from  the  camera,  (j>j^(r),  and  the  simulated  photon  flux,  0^(r)  at  identical  sample  locations.  This  is 
accomplished  by  minimizing  the  L2-norm  of  the  objective  function  in  an  optimization  routine: 

^  “  {^s  ~ 

where  Q  is  the  objective  function  to  minimize.  This  is  an  underdetermined  problem,  as  measurements 
are  made  only  at  the  boundary.  Because  this  problem  is  underdetermined,  the  model,  G,  is  linearized 
with  a  Taylor  approximation  and  formed  into  an  iterative  algorithm  as: 

G(c-)  =  G(Cj_y)  +  G' A.C 

where  G’  is  the  partial  differential  of  the  model  with  respect  to  the  concentration,  also  known  as  the 
Jacobian,  J.  Minimizing  Q  with  respect  to  c  and  substituting  GlcJ  into  and  J  for  G’  yields: 
27{(G(c,_i)  +  7Ac)-(/)^)  =  0 
Solving  for  the  concentration  yields: 

This  problem  is  ill-posed,  so  it  is  solved  using  the  Levenberg-Marquardt  algorithm,  which  includes  a 
stabilization  parameter.  A,  in  the  inversion  to  avoid  singularities: 

Ac  =  -[yV+A;]''y"(0,-0„) 

This  is  iterated  until  a  minima  is  reached  (the  L2  norm  of  the  update  is  less  than  1%  of  the  previous 
iteration),  or  until  15  iterations  are  performed,  whichever  occurs  earlier.  The  stabilization  parameter  is 
reduced  at  each  iteration  as  the  algorithm  approaches  the  minimum  and  converges  on  the  solution. 

The  hybrid  X-ray  /  diffuse  optical  model  was  demonstrated  in  a  breast-sized  phantom,  simulating 
a  breast  lumpectomy  geometry.  Both  numerical  and  experimental  phantoms  were  tested,  with  lesion- 
simulating  objects  of  various  sizes  and  depths.  Results  showed  localization  accuracy  with  median  error 
of  2.2mm,  or  4%  of  object  depth,  for  small  2-14mm-diameter  lesions  positioned  from  1cm  to  4.5cm  in 
depth.  This  compares  favorably  with  fluorescence  optical  imaging,  which  is  not  able  to  resolve  such 
small  objects  at  this  depth.  The  recovered  lesion  size  had  lower  size-bias  in  the  X-ray  excitation 
direction  than  the  optical  direction,  which  was  expected  due  to  the  optical  scatter.  However,  shown  in 
Figure  4,  the  technique  was  shown  to  be  quite  invariant  in  recovered  size  with  respect  to  depth,  as  the 
standard  deviation  was  less  than  2.5mm.  Sensitivity  was  a  function  of  dose  (shown  in  Figure  4a); 
radiological  doses  were  found  to  provide  sufficient  recovery  for  pg/ml  concentrations,  while  therapy 
dosages  provided  recovery  for  ng/ml  concentrations.  Experimental  phantom  results  agreed  closely  with 
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the  numerical  results,  with  positional  errors  recovered  within  8.6%  of  the  effective  depth  for  a  5mm 
object,  and  within  5.2%  of  the  depth  for  a  10mm  object.  Object  size  median  error  was  within  2.3%  and 
2%  for  the  5mm  and  10mm  objects,  respectively.  For  shallow-to-medium  depth  applications  where 
optical  and  radio-emission  imaging  modalities  are  not  ideal,  such  as  in  intra-operative  procedures,  this 
new  technique,  LAXLT,  may  be  a  useful  tool  to  detect  molecular  signatures  of  disease.  For  this  aim, 
the  software  and  hardware  proof-of-principle  experiments  were  performed. 


a 


H —  1  pg/ml 
■* —  1  ng/ml 
« —  1  ng/ml 
1  mg/ml 


b 


-  e  -  Truth:  5:1 
— e —  Expt:  5:1 

-  B  -Truth:  10:1 
^^Expt:  10:1 

-  0  -  Truth:  25:1 

Expt:  25:1 

-  V  -  Truth:  40:1 
V  Expt:  40:1 


Object  Depth  (mm) 


Figure  4:  (a)  Numerical  phantom  experiment  examining  the  effect  of:  (a)  variable  object  to  background 
contrast  and  (b)  variable  concentration,  of  a  6mm-diameter  object  vs.  depth  from  the  detection  plane. 


2.3:  Imaging  Hardware  Training 

In  performing  these  tasks,  the  investigator  explored  various  methods  of  multimodality  imaging,  by 
fusing  X-ray  and  optical  photon  propagation  models.  This  work  was  carried  out  under  guidance  from 
Dr.  Lei  Xing,  as  proposed,  who  is  an  expert  in  X-ray  modeling  in  tissue.  A  numerical  photon 
propagation  model  was  built,  and  integrated  into  the  PFs  existing  code  for  optical  photon  modeling. 
The  system  development  for  pre-clinical  imaging  was  a  first  for  this  PI:  a  small-animal  system  devoted 
to  testing  modality  feasibility,  and  one  that  may  have  potential  for  high-resolution  small-animal 
imaging. 


2.4:  Imaging  Software  Multiplexing  Development 

The  ability  to  image  multiple  nanoparticle-based  contrast  agents  simultaneously  is  a  key  advantage  of 
RLI  over  other  molecular  imaging  modalities,  such  as  PET,  and  many  types  of  fluorescence  imaging. 
Since  the  radioluminescent  nanoparticles  emit  optical  light  at  unique  wavelengths,  depending  on  their 
lanthanide  dopant,  we  are  able  to  selectively  image  the  unique  phosphors  with  the  use  of  optical  filters. 
We  demonstrated  this  technique  by  separating  two  distinct  nanophosphor  contrast  agents  in  gelatin 
phantoms  and  a  small  animal  phantom. 

To  enable  multiplexed  imaging,  both  hardware  and  software  were  developed.  A  filter  wheel 
(Caliper  Life  Sciences,  Alameda,  CA)  was  mounted  on  our  imager  with  2  filters  for  2  different  phosphor 
nanoparticles  (548nm:  CaF2:Tb,  700nm:  CaF2:  Eu).  Leaded  glass  was  placed  above  the  subject.  To 
form  images  of  the  concentrations,  c,  of  each  of  these  particles,  the  light  emitted,  (j),  was  recorded  over 
the  several  optical  spectral  regions  defined  by  filter  fi  and  filter  f2  to  form  images  of  the  subject,  ij  =^fi) 
and  i2  =(p(fj)-  These  images,  and  the  pre-recorded  reference  spectra,  e,  for  Bao.55Yo.3E2:Tb^'^,  (erb). 
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Bao.ssYo.aFiiEu^'*’,  {seu),  were  input  into  a  linear-least  squares  algorithm  to  extract  the  contributions  from 
each  nanophosphor.  Reference  spectra  were  collected  with  a  calibrated  spectrometer,  thus  enabling  a 
quantitative  comparison  between  the  nanophosphors. 

To  validate  the  ability  of  this  method  to  image  multiple  nanophosphors  simultaneously,  five 
phantoms  (99%  DI  water,  1%  agarose)  were  fabricated  containing  Tb^"^  and  Eu^'*'  doped  RENPs  of 
linearly  varying  concentrations  (in  units  of  mg/ml)  and  mixed  in  the  following  ratios:  (0:10;  2. 5:7. 5;  5:5; 
7. 5:2. 5;  and  10:0).  Eigure  5(a)  displays  the  unmixed  images  for  linearly  increasing  (top-to-bottom) 
amounts  of  Tb^^-doped  RENPs  on  the  left,  and  decreasing  (top-to-bottom)  Eu^'^-doped  RENPs  on  the 
right,  at  each  concentration.  The  median  values  for  each  region  of  interest  were  plotted  with  respect  to 
concentration.  Eigure  5(b)  shows  the  median  raw  signal  detected  from  each  nanophosphor 
concentration.  Eigure  5(c)  demonstrates  this  method’s  efficacy  in  separating  the  nanophosphors;  the 
correlation  between  the  samples  was  (r  =  -0.98). 
a  BaYF4:Tb 
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Figure  5:  (a)  Shown  are  the  Bao.55Yo.3E2:Tb^'^  and  Bao.55Yo.3E2:Eu^'^  phosphor  concentrations  with 
increasing/decreasing  concentration  (top-to-bottom)  of  Bao.55Yo.3E2:Tb^'^  /  Bao.55Yo.3E2:Eu^'^, 
respectively,  (b)  The  raw  signal  detected  for  each  respective  nanophosphor,  (c)  The  relative  median 
recovered  concentration  in  each  ROI  plotted  with  respect  to  concentration. 

To  demonstrate  this  method  in  a  pre-clinical  mouse  model,  four  batches  of  RENPs  were  each  mixed 
with  52  pCi  (Curie)  of  activity  of  18-P  and  Matrigel  (BD  Biosciences,  Sparks,  MD)  and  injected  into 
both  forelegs  and  flanks  of  a  euthanized  nude  mouse,  with  Eu^'^-doped  RENPs  (Eu)  only  on  the  left 
foreleg,  Tb^’^-doped  RENPs  (Tb)  only  on  the  left  flank,  an  equal  mixture  of  Eu^'^-doped  RENPs  and 
Tb^'^-doped  RENPs  (Mix)  on  the  right  flank,  and  an  inactive  undoped  nanophosphor  control  (Ctrl)  on  the 
right  foreleg.  The  locations  of  the  subcutaneously  injected  RENPs  are  shown  in  Eigure  6(a).  The 
feasibility  of  this  method  in  a  pre-clinical  subject  with  X-ray  excitation  is  shown  in  Eigure  6(b)-(d). 
Tb^'^-doped  RENPs  embedded  in  the  left  and  right  flanks  are  shown  successfully  recovered  in  Eigure 
6(b).  The  Eu^'^-doped  RENPs  in  the  left  foreleg  and  the  right  flank  are  shown  successfully  recovered  in 
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Figure  6(c).  As  a  control,  inactive  RLNPs  in  the  right  foreleg  did  not  luminesce.  Figure  6(d)  shows  the 
results  of  the  spectral  unmixing  algorithm,  with  the  Eu^’^-doped  RLNPs  shown  in  red,  and  the  Tb^"^- 
doped  RLNP  shown  in  green.  The  yellow  in  the  multiplexed  image  indicates  the  presence  of  both  types 
of  nanophosphors. 


Figure  6:  X-ray  Luminescence:  (a)  Schematic  of  the  locations  of  each  type  of  RLNP.  The  inactive 
particles  are  indicated  with  the  abbreviation  (Ctrl).  (b,c)  The  unmixed  signal  from  the  Bao.55Yo.3p2:Tb 
and  Bao.55Yo.3p2:Eu  particles,  respectively,  (d)  The  unmixed  multiplexed  image  with  colorbars  for  the 
relative  concentrations  of  Bao.55Yo.3p2:Tb  and  Bao.55Yo.3p2:Eu. 


2.5;  Imaging  Hardware:  Biopsy  Needle  -  based  RLI 

Currently,  biopsy  is  guided  by  X-ray  fluoroscopy  or  ultrasound;  these  imaging  modalities  do  not  offer 
the  ability  to  image  molecular  contrasts  of  tissue,  and  thus  are  limited  in  their  ability  to  localize  cancer. 
This  shortcoming  is  especially  true  in  small  cancerous  lesions,  where  structure  visible  by  X-rays  or 
ultrasound  is  not  yet  established;  there  is  no  current  suitable  means  to  localize  cancers  at  this  critical 
stage. 

To  meet  this  need,  as  outlined  in  Aim  2,  a  custom  biopsy  RLI  catheter  was  fabricated  to 
demonstrate  proof-of-principle  for  phosphor  image-guided  biopsy  procedures.  The  specifications  were 
to  provide  a  <  10  gauge  imaging  instrument  that  would  fit  down  the  needle  bore.  The  ability  to  image  is 
paramount  to  enable  a  visualization  of  the  target  mass,  and  to  enable  the  search  for  the  molecular 
beacon.  The  schematic  of  this  instrument,  shown  in  Ligure  7,  utilizes  a  water-cooled  electron  multiplied 
CCD  (ProEM,  Princeton  Instruments)  attached  to  the  appropriate  4: 1  magnification  lenses  (35mm  and 
8mm,  Edmunds  Optics,  Megapixel)  to  project  the  proximal  end  of  a  leached  fiber  bundle  (Schott,  Inc.) 
onto  the  full  CCD  field  of  view.  This  fiber  bundle  is  protected  by  a  heat  shrink  catheter  made  of  PTLE 
to  prevent  kinking  and  breaking  of  the  individual  fibers. 

Images  of  the  detection  hardware  and  lens-coupled  fiber  are  shown  in  Ligure  8.  The  biopsy 
imaging  system  was  tested  for  focus  by  imaging  a  business  card  under  white-light  illumination;  an 
example  image  is  shown  in  Ligure  8c. 
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Figure  7:  Schematic  of  the  RLI  biopsy  device,  consisting  of  a  cooled  CCD,  a  filter  wheel  for 
multiplexing  ability,  and  a  lens  at  the  distal  tip. 


Figure  8:  (a,b)  Photographs 
of  the  Radioluminescent 
Biopsy  Imaging  device, 
specifically  of  (a)  the  CCD 
camera,  filter  wheel,  and 
coupling  optics;  (b)  a  close-up 
of  the  distal  tip  of  the  catheter, 
(c)  white-light  image  of  text 
from  a  business  card  taken 
with  the  RLI  biopsy  device. 


2.6:  Imaging  Hardware:  RLI  Endoscopy 

With  a  similar  goal  for  interventional  imaging,  a  larger  bore  fiber-based  system  was  built  utilizing  an 
imaging  optical  fiber  coupled  to  a  highly  sensitive  ICCD  camera.  This  system  had  a  larger  aperture,  and 
thus  was  more  sensitive  to  Radioluminescence.  Images  of  the  system  are  shown  in  Figure  9  a-b.  A 
cartoon  of  the  application  of  this  system  for  surgical  oncology  guidance  is  detailed  in  Figure  9c. 
Specifically,  a  micro-imaging  lens  (Schneider  Cinegon,  F/1.4,  12mm  focal  length)  was  coupled  at  the 
distal  end  of  an  optical  imaging  fiber  bundle  that  was  108  mm  long,  with  a  5  x  6.7  mm  active  area.  The 
pixels  of  the  fiber  bundle  are  made  of  discrete  10pm  fibers  (Schott,  IG154).  At  the  proximal  end,  a 
Pentax  (F/1.4)  lens  provided  relay  optics  to  the  camera.  The  camera,  an  image-intensified  CCD 
(Stanford  Photonics  Inc.,  Turbo  640-Z),  had  640  x  480  pixels  and  single-photon  imaging  capability.  To 
minimize  background  light  and  to  emulate  a  light-tight  anatomical  cavity,  all  images  were  taken  in  a 
dark  box.  In  addition,  images  were  post-processed  with  a  software  thresholding  technique  that  removed 
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pixel  values  below  a  manufaeturer-reeommended  value;  this  optimal  value  was  ehosen  to  remove  low- 
intensity  pixels  whieh  corresponded  to  thermal  and  read  noise  on  the  CCD.  These  steps  reduced  the 
noise  to  -20  counts  per  second. 


Figure  9.  A  prototype  of  fiber  based  system  for  endoscopic  and  laparoscopic  Cerenkov  Luminescence 
Imaging  applications.  An  optical  fiber  with  micro-imaging  lens  (a)  was  coupled  to  an  image-intensified 
CCD  (b).  (c)  A  suggested  application  for  the  fiber  based  system  for  endoscopic  and  laparoscopic 
Radioluminescence  Imaging. 

We  investigated  the  feasibility  of  RNLP-free  Radioluminescence  Imaging  with  the  interventional  system 
shown  in  Figure  10  for  guiding  cancer  imaging  by  performing  spatial  resolution,  sensitivity,  and  proof- 
of-concept  experiments. 

The  spatial  resolution  of  the  system  was  characterized  using  a  standard  PET/SPECT  phantom 
(Micro  Hot-spot  phantom.  Data  Spectrum  Corporation,  Hillsborough,  NC)  filled  with  15.2  MBq  (410 

1  o 

pCi)  of  E-EDG;  this  phantom  is  also  commonly  used  to  evaluate  the  spatial  resolution  of  PET  systems. 
Both  the  2.4  and  1.6  mm  cylindrical  holes  were  used  to  form  line  profiles,  taken  from  the  fiber-based 
images.  Two  images  were  acquired:  an  ambient-light  image,  and  a  functional  CEI  image.  The  exposure 
time  for  all  ambient  images  was  <1  sec,  while  the  functional  image  acquisition  time  was  5  minutes.  The 
subject  was  located  5cm  from  the  optical  system.  The  photographic  and  Cerenkov  images  taken  using 
the  standard  optical  imaging  system  are  similar  to  the  photographic  and  Cerenkov  images  taken  with  the 
fiber-based  system  (Eigure  lOA).  The  line  profiles  demonstrate  a  high  correlation  between  the  Cerenkov 
and  photographic  line  profiles.  Eor  the  2.4  mm  holes  (Eigure  lOB),  the  line  profile  limits  depicted  by  the 
large  green  stars  in  Eigure  lOA,  had  a  Pearson  correlation  coefficient  of  r  =  0.78  (significance:  p  <  le-5). 
Eor  the  1.6  mm  holes,  the  line  profile  limits  depicted  by  the  smaller  red  stars,  the  correlation  between  the 
line  profiles  of  r  =  0.71  (significance:  p  <  le-5).  The  peaks  were  also  in  near  identical  locations;  the 
difference  in  peak  locations  between  the  holes  for  the  line  profiles  of  the  photographic  and  Cerenkov- 
emission  images  was  3.8%,  a  difference  of  only  a  few  pixels. 

Eigure  lOE  shows  the  reduction  in  signal  during  E-EDG  decay  over  time  for  two  wells:  one 
filled  with  E-EDG,  and  one  filled  with  a  water/glycerol  mixture  only.  With  a  5  min  integration  time,  a 
minimum  of  -45  KBq  (1.21  pCi)  of  activity  can  be  identified  as  different  than  the  control  well 
containing  the  water/glycerol  solution  (signal  to  noise  ratio  >  1).  Note  the  decrease  in  signal  in  the 
control  well;  this  change  in  signal  was  due  to  stray  gamma  photons  from  the  E-EDG  well  interacting 
with  the  optical  fiber  and  creating  scintillation  background  light  in  the  image. 
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Figure  10.  Characterization  of  the  fiberseopie 
system  spatial  resolution.  (A)  Respeetive  ambient 
(left)  and  Radiolumineseenee  images  (right)  of  a 
(4.3  cm  inner  diameter  phantom).  (B-E) 
Quantitative  line  profiles  of  the  ambient  and 
Cerenkov  images  for  the  1.6mm  &  2.4mm  holes, 
respectively;  the  line  sampled  is  indicated  by  the 
boxes  in  (A).  (F)  Sensitivity  with  FDG. 


2.7;  Imaging  Hardware  Development 

In  performing  these  tasks,  the  investigator  gained  valuable  training  in  optical  systems  and  catheter-based 
devices.  This  experience  will  be  useful  for  future  device  development.  In  addition,  the  investigator 
gained  valuable  experienee  in  working  with  small-animal  phantoms  (post-mortum)  for  more  biologieally 
appropriate  phantom  demonstration,  where  the  tissue  properties  and  geometries  mateh  that  of  the  living 
subjects. 


3.1:  Imaging  Hardware  Development:  Sensitivity  Comparison 

Towards  meeting  aim  3,  we  investigated  the  sensitivity  of  RLI  and  compared  it  to  the  current  state-of- 
the-art  in  fluorescenee  molecular  imaging.  To  capture  sensitivity  for  the  RLI  system,  we  imaged 
decreasing  concentrations  of  CaF2:Eu,  coated  with  Au,  phosphor  mixed  with  1%  agar.  As  shown  in 
Figure  I  la,  we  were  able  to  linearly  reeover  eoneentrations  of  phosphor  to  lOpg/ml,  or  about  15nM. 

This  eompares  favourably  to  fluoreseenee  imaging,  shown  in  Figure  11b  [1],  whieh  has  a  sensitivity 
limitation  of  about  5nM. 

However,  RLI  ean  be  improved.  As  is  apparent  in  Figure  I  la,  our  investigations  have  uneovered 
a  potential  limit  to  sensitivity.  While  irradiating  water-soluble  phosphor  samples,  we  discovered  that  X- 
rays  create  optical  luminescence  in  water  [2,3],  which  has  been  attributed  to  hydroxide  radicals  and 
hydrated  electrons  in  the  literature.  This  signal  adds  baekground  signal  to  RLI,  potentially  limiting 
sensitivity  to  20pg/ml,  or  30nM,  as  shown  below. 
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Figure  11:  (a)  Sensitivity 
of  Radioluminescence 
imaging.  The  signal  from 
several  phantoms  with 
concentrations  ranging 
from  lOpg/ml  to  10  mg/ml 
is  shown  in  black.  The 
background  water 
luminescence  is  shown  in 
red.  (b)  Sensitivity  of 
fluorescence  imaging  for 
two  different  filtering 
strategies  (from  [1]).  The 
fluorescence  sensitivity  is 
limited  by  the 
autofluorescence  in  the 
tissue,  or  in  this  case,  the 
autofluorescence  of  lipid  in 
the  intralipid  /  Indocyanine 
green  phantom. 


4.1:  Feasibility  of  injectable  nano-scintillators  in  vivo 

Significant  effort  was  dedicated  to  RLNP  Radioluminescence  Imaging;  work  is  ongoing.  Due  to  the 
difficulty  in  fabricating  biocompatible,  stable  nanoparticles,  progress  was  hindered.  In  addition,  during 
our  investigation  of  illumination  sources,  we  noticed  that  the  light  emitted  from  injected  FDG  was 
nearly  as  bright  as  that  from  RLI  nanoparticles  (see  Figure  12),  depending  on  the  efficiency  of  the 
nanophosphor.  After  literature  investigation,  it  emerged  that  the  light  emitted  from  the  radiotracer  was 
Cerenkov  light.  As  shown  in  Figure  12,  the  Cerenkov  light  was  able  to  be  unmixed  from  the  RLNP 
signal.  Thus,  an  investigation  into  the  utility  of  Radioluminescence  (Cerenkov)  imaging  for 
interventional  imaging  was  launched,  utilizing  the  same  instrumentation  designed  for  RLNP  imaging. 
We  believe  that  RLNP-free  RLI  imaging  has  significant  potential  in  several  niche  clinical  applications, 
such  as  surgical  guidance  for  the  removal  of  breast  cancer.  Below  are  our  investigations  into  RLNP-free 
RLI. 
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Figure  12:  X-ray  Luminescence: 
(a)  Schematic  of  the  locations  of 
each  type  of  RLNP.  The  inactive 
particles  are  indicated  with  the 
abbreviation  (Ctrl).  (b,c)  The 
unmixed  signal  from  the 
Bao.55Yo.3F2:Tb  and 

Bao.55Yo.3F2:Eu  particles, 

respectively,  (d)  The  unmixed 
multiplexed  image  with  colorbars 
for  the  relative  concentrations  of 
Bao.55Yo.3F2:Tb  and 

Bao.55Yo.3F2:Eu. 


4.2:  Feasibility  of  RLNP-free  Radioluminescence  imaging 

For  this  investigation,  we  studied  the  feasibility  of  RNLP-free  Radioluminescence  Imaging  with  the 

1  Q 

interventional  system  shown  in  Figure  13  for  guiding  in-vivo  cancer  imaging.  Specifically,  F-FDG 
was  produced  by  Radiochemistry  Facility  at  Stanford  University  (Stanford,  CA).  Rat  glioma  cell  line  C6 
was  obtained  from  American  Type  Culture  Collection  (Manassas,  VA).  Female  athymic  nude  mice 
(nu/nu)  (N=5)  obtained  from  Charles  River  Laboratories,  Inc.  (Cambridge,  MA)  were  4-6  weeks  of  age. 
A  CRC-15R  PET  dose  calibrator  (Capintec  Inc.,  Ramsey,  NJ)  was  used  for  all  radioactivity 
measurements.  All  animal  studies  were  carried  out  in  compliance  with  federal  and  local  institutional 
guidelines  for  the  conduct  of  animal  experimentation.  C6  cells  were  cultured  in  DMEM  medium 
supplemented  with  10%  fetal  bovine  serum  (EBS)  and  1%  penicillin/streptomycin  (Invitrogen  Eife 
Technologies,  Carlsbad,  CA).  The  cell  line  was  maintained  in  a  humidified  atmosphere  of  5%  CO2  at 
37°C,  with  the  medium  changed  every  other  day.  A  75%  confluent  monolayer  was  detached  with  trypsin 
and  dissociated  into  a  single  cell  suspension  for  further  cell  culture.  Approximately  1  x  10^  C6  cells 
suspended  in  phosphate  buffered  saline  (PBS,  O.IM,  pH=7.2,  Invitrogen,  Carlsbad,  CA)  were  implanted 

a 

subcutaneously  in  the  left  legs  of  nude  mice.  Tumors  were  allowed  to  grow  to  a  size  of  150  to  200  mm 
(2-3  weeks),  and  the  tumor-bearing  mice  were  imaged  in  vivo.  Validation  was  performed  with  an  IVIS 
Spectrum  system  (Caliper  Eife  Science,  Hopkinton,  MA).  Eor  all  in  vivo  studies,  radionuclides  were 
diluted  in  PBS.  Animals  were  placed  in  a  light-tight  chamber  under  isofluorane  anesthesia.  Each 
acquisition  took  3  min  for  all  studies  without  filters.  Images  were  acquired  and  analyzed  using  Eiving 

1  Q 

Image  3.0  software  (Caliper  life  sciences,  Hopkinton,  MA).  The  mice  were  fasted  overnight  prior  to  E- 
EDG  imaging  and  kept  anesthetized  by  inhalation  of  2%  isofluorane  during  the  experiment. 

The  REI  images  comparing  both  systems  before  removal  of  the  tumor  for  mouse  1  are  shown  in 
Eigure  13;  the  IVIS  system  produced  Eigures  1 3 A,  while  the  fiber-based  CEE  system  produced  Eigures 
13B.  Comparing  images  in  Eigure  13B  quantitatively,  the  Cerenkov  signal  in  the  tumor  to  background 
for  the  excised  tumor  tissue  for  the  fiber-based  CEE  system  was  determined  by  computing  a  ratio 
between  the  median  value  in  the  region  of  interest  (ROI)  encircling  the  tumor  to  the  median  value  in  the 
cleared  tumor  cavity.  Eor  the  removed  tumor,  the  tumor-to-background  was  1.28  (for  comparison,  the 
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tumor-to-background  in  the  IVIS  was  1.16  for  the  3  minute  sean).  For  Mouse  2-5,  the  tumor-to- 
background  for  the  removed  tumor  was  1.41,  1.21,  1.02,  and  1.17,  respeetively.  Tumor  tissue  light 
emission  was  significantly  higher  than  the  exposed  cavity  for  all  mice  (Student  t-test  for  paired  samples, 
p  <  0.05  for  all).  Much  of  the  low  tumor-to-background  value  in  Mouse  4  can  be  explained  by  the 
reflection  of  the  light  emitted  from  the  tumor  by  the  tumor  cavity  which  was  directly  adjacent.  Residual 
tumor  tissue  after  surgery  is  also  possible. 


Figure  13.  Mouse  1  imaged  by  (A)  the  IVIS  optical  system,  consisting  of  an  ambient- light  image  (left), 
the  luminescent  image  (middle),  and  the  fused  image  (right);  and  (B)  fiber  based  system  after  surgery  to 
remove  tumor  tissues,  consisting  of  an  ambient-light  image  (left),  the  luminescent  image  (middle),  and 
the  fused  image  (right).  (C)  Comparison  of  tumor  to  background  for  all  mice  in  the  study. 


4.3:  Small-Animal  in-vivo  Imaging  Training 

In  performing  these  tasks,  the  investigator  gained  valuable  training  in  handling  small-animal  mouse 
models  and  performing  in-vivo  experiments.  This  experience  will  be  useful  for  future  molecular 
imaging  ventures. 


KEY  RESEARCH  ACCOMPLISHMENTS 


•  Established  the  feasibility  of  Radioluminescent  Tomography  in  reflection  mode. 

•  Developed  a  hybrid  X-ray /Optical  photon  propagation  model  to  perform  high  resolution 
Radioluminescent  Tomography,  in  an  appropriate  geometry  for  Intraoperative  Radiation 
Therapy. 

•  Characterized  nanophosphors  suitable  for  RET,  including  Gadolinium  oxy sulfide  and  Barium 
Yttrium  Eluoride. 

•  Eabricated  a  small-animal  pre-clinical  RET  imaging  box  to  enable  automated,  controlled,  RET 
experiments. 

•  Improved  instrumentation  through  leaded  glass  and  emission  filters  to  enable  sensitivity  of 
Radioluminescent  Imaging  to  30  nM. 
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•  Developed  strategy  to  remove  water- luminescent  component  to  improve  sensitivity. 

•  Demonstrated  multiplexing  ability  of  RLI  in  phantoms  and  small  animal  phantoms,  and 
completed  both  hardware  and  software  requirements. 

•  Fabricated  two  imaging  fiber-based  RLI  systems,  including  a  biopsy  system  and  a  hand-held 
system  to  demonstrate  proof-of-concept. 

•  Fully  characterized  an  RLNP-free  RLI  (Cerenkov)  imaging  system  in  terms  of  sensitivity  and 
spatial  resolution. 

•  Demonstrated  RLNP-free  RLI  (Cerenkov)  imaging  in  successfully  removing  cancerous  tissue 
from  a  small- animal  mode. 

REPORTABLE  OUTCOMES; 

12  journal  publications;  20  conference  abstracts,  and  9  conference  presentations;  several  postdoctoral 
fellows  are  now  being  mentored  and  working  in  RLI;  4  courses  taken,  including  BioE222:  Molecular 
Imaging,  Med374:  Medical  Device  Design  (taken  after  the  grant  was  accepted,  before  it  was  funded), 
the  Comprehensive  Cancer  Training  Program,  and  SIE:  The  Stanford  Institute  for  Entrepreneurs,  and  3 
patent  applications. 
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CONCLUSION 

This  grant  has  enabled  the  grantee  to  focus  on  the  task  of  developing,  characterizing,  and  implementing 
new  molecular  imaging  strategies,  and  perform  imaging  in  vivo.  In  particular,  the  Radioluminescence 
Imaging  modality  was  thoroughly  investigated.  The  incorporation  of  hardware  and  software  innovations 
has  resulted  in  high  sensitivity  which  is  demonstrated  as  suitable  for  small  animal  imaging.  The 
additional  innovation  of  multiplexed  imaging  methods  will  aid  the  research  and  development  of  RLI  of 
RLNPs,  by  enabling  multiple  molecular  targets  to  be  simultaneously  imaged.  In  addition,  this  funding 
period  has  resulted  in  experience  in  mentoring  other  fellows. 

The  research  and  training  completed  are  significant  for  the  eradication  of  breast  cancer,  and  for 
training  for  future  novel  tools  for  early  detection  of  breast  cancer.  The  coursework  and  experience 


21 


gained  by  this  PI  will  be  critical  in  his  goal  to  translate  breast  cancer  treatment  technology  into  the 
clinic. 
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Abstract: 

X-ray  Luminescence  Tomography  (XLT)  has  recently  been  proposed  as  a  new  imaging 
modality  for  biological  imaging  applications.  This  modality  utilizes  phosphor 
nanoparticles  which  luminesce  near-infrared  light  when  excited  by  X-ray  photons.  The 
advantages  of  this  modality  are  that  it  uniquely  combines  high  sensitivity  of 
radioluminescent  nanoparticles  and  high  spatial  localization  of  collimated  X-ray  beams. 
Currently,  XLT  has  been  demonstrated  using  X-ray  spatial  encoding  to  resolve  the 
imaging  volume.  However,  there  are  applications  where  the  X-ray  excitation  may  be 
limited  by  geometry,  where  increased  temporal  resolution  is  desired,  or  where  a  lower 
dose  is  mandatory.  This  paper  extends  the  utility  of  XLT  to  meet  these  requirements  by 
incorporating  a  photon  propagation  model  into  the  reconstruction  algorithm  to  recover 
dimensions  remaining  in  an  X-ray  limited  angle  (LA)  geometry.  This  enables  such 
applications  as  image-guided  surgery,  where  the  ability  to  resolve  lesions  at  depths  of 
several  centimeters  can  be  the  key  to  successful  resection.  The  hybrid  X-ray  /  diffuse 
optical  model  is  first  formulated  and  then  demonstrated  in  a  breast-sized  phantom, 
simulating  a  breast  lumpectomy  geometry.  Both  numerical  and  experimental  phantoms 
are  tested,  with  lesion-simulating  objects  of  various  sizes  and  depths.  Results  show 
localization  accuracy  with  median  error  of  2.2mm,  or  4%  of  object  depth,  for  small  2- 
14mm-diameter  lesions  positioned  from  1cm  to  4.5cm  in  depth.  This  compares  favorably 
with  fluorescence  optical  imaging,  which  is  not  able  to  resolve  such  small  objects  at  this 
depth.  The  recovered  lesion  size  has  lower  size-bias  in  the  X-ray  excitation  direction 
than  the  optical  direction,  which  is  expected  due  to  the  increased  optical  scatter. 

However,  the  technique  is  shown  to  be  quite  invariant  in  recovered  size  with  respect  to 
depth,  as  the  standard  deviation  is  less  than  2.5mm.  Sensitivity  is  a  function  of  dose; 
radiological  doses  are  found  to  provide  sufficient  recovery  for  pg/ml  concentrations, 
while  therapy  dosages  provide  recovery  for  ng/ml  concentrations.  Experimental  phantom 
results  agree  closely  with  the  numerical  results,  with  positional  errors  recovered  within 
8.6%  of  the  effective  depth  for  a  5mm  object,  and  within  5.2%  of  the  depth  for  a  10mm 
object.  Object  size  median  error  is  within  2.3%  and  2%  for  the  5mm  and  10mm  objects, 
respectively.  For  shallow-to-medium  depth  applications  where  optical  and  radio- 


emission  imaging  modalities  are  not  ideal,  such  as  in  intra-operative  procedures,  LAXLT 
may  be  a  useful  tool  to  detect  molecular  signatures  of  disease. 

1.  Introduction 

Imaging  plays  a  vital  role  in  the  management  of  cancer  care,  for  detection,  staging, 
intervention,  and  monitoring  of  treatment  response.  Despite  its  ubiquitous  use  elsewhere, 
the  role  of  imaging  in  surgery  is  limited,  as  it  is  dominated  by  C-arm  fluoroscopy  and 
optical  endoscopy.  These  tools  are  appropriate  for  visualizing  tissue  structure,  yet  are 
limited  in  their  sensitivity  to  microscopic  disease.  This  limitation  affects  such  procedures 
as  surgical  breast  lumpectomy,  as  many  studies  have  found  that  surgeons  are  unable  to 
remove  all  tumor  tissue  present  in  the  surgical  field  (for  example,  Gibson  et  al.  identified 
residual  tumor  in  55%  of  the  cases[l]).  The  risks  of  local  failure  are  high,  as  local  failure 
often  leads  to  distant  metastasis  [2].  Thus,  there  is  a  need  for  tools  to  provide  surgeons 
with  more  sensitive,  more  specific  image-guidance. 

This  need  may  be  fulfilled  with  molecular  imaging,  which  promises  to  image 
molecular  and  cellular  processes,  and  may  allow  the  early  identification  of  disease  or 
status  of  disease  progression  and  treatment  [3].  Developing  these  tools  for  the  operating 
room  would  aid  a  physician  during  an  intervention,  by  allowing  the  clinician  to  identify 
near-microscopic  regions  of  disease,  such  as  at  the  tumor  margin.  Ideally,  this  tool  would 
be  able  to  image  at  a  depth  of  several  centimeters,  so  that  disease  buried  beneath  the 
superficial  layers  could  be  identified.  Several  potential  applications  for  this  technology 
could  be  in  removing  occult  disease  in  breast[4,  5],  brain[6]  and  hepatic  tumors[7],  where 
imaging  is  currently  being  incorporated  into  the  clinic,  and  new  innovations  may  be 
readily  translated. 

This  paper  develops  and  demonstrates  a  novel  X-ray  Luminescence  Tomographic 
(XLT)  method  that  is  uniquely  suited  for  image-guided  surgical  applications.  This 
method.  Limited- Angle  XLT  (LAXLT),  utilizes  a  photon  propagation  model  to  enable 
XLT  for  surgical  guidance,  where  XLT’s  advantages  are  the  clearest  for  translation  into 
the  clinic.  XLT  has  been  recently  introduced[8]  and  demonstrated  in  simulation  and  in 
phantoms  [9,  10].  This  imaging  modality  utilizes  nano-sized  phosphors  which  emit 
optical  near-infrared  light  upon  X-ray  excitation  [11,  12].  Attaching  these  phosphors  to 


molecular  probes  (e.g.  antibodies  and  peptides)  that  target  molecular  markers  specific  to 
tumors,  such  as  angiogenesis  markers  like  epidermal  growth  factor  receptor [13],  or  a^Pj- 
integrin  [14,  15]  expression  could  allow  the  surgeon  to  differentiate  between  normal  and 
cancerous  tissue.  XLT  has  several  advantages  to  current  molecular-sensitive  imaging 
modalities:  emission  imaging  techniques,  such  as  gamma  cameras,  are  limited  in  their 
ability  to  discriminate  depth  due  to  the  limited  angles  that  may  be  imaged  during  surgical 
procedures [16];  optical  imaging,  on  the  other  hand,  has  the  ability  to  provide  depth 
localization,  and  is  currently  under  investigation  for  surgical  guidance  [6,  4,  17],  yet  is 
limited  in  its  ability  to  image  deeper  than  ~lcm[18].  Depth  is  important  to  discern  to 
determine  occult  lesions  lying  under  the  superficial  layer,  and  to  determine  the  feasibility 
of  surgical  removal  of  a  lesion. 

XLT  utilizes  the  extremely  low-scatter  of  X-rays  compared  to  optical  fluorescence 
imaging  to  enable  higher  spatial  resolution.  A  thin  pencil-beam  of  collimated  X-rays 
may  be  maintained  while  the  X-rays  propagate  through  tissue  of  several  cm;  this  spatial 
localization  is  in  contrast  to  optical  excitation,  which  is  highly  attenuated  and 
scattered[19].  By  rotating  the  X-ray  (or  similarly  the  phantom)  to  cover  all  angular 
projections,  the  resolving  power  is  limited  merely  by  the  width  of  the  beam  (up  to  the 
diffraction  of  the  X-ray).  A  numerical  analysis  demonstrated  that  2.25cm  deep  objects  as 
small  as  1mm  (using  a  1mm  beam  width)  with  a  nanoparticle  concentration  of  0.4pM 
could  be  resolved;  incraesing  dose  increased  the  sensitivity [9].  However,  there  are 
applications  where  the  X-ray  excitation  may  be  limited  by  geometry,  where  increased 
temporal  resolution  is  desired,  or  where  a  lower  dose  is  mandatory;  one  such  application 
is  intraoperative  breast  cancer  lumpectomy,  where  it  may  not  be  possible  or  desirable  to 
irradiate  over  the  full  projection  space.  In  these  cases,  it  would  be  beneficial  to  irradiate 
over  a  limited  projection  space,  and  use  the  ability  of  the  optical  detectors  to  resolve  the 
remaining  dimensions.  This  technique  could  also  have  utility  in  decreasing  dose  to  the 
tissue,  as  fewer  irradiation  beamlets  are  needed  to  resolve  the  volume. 

This  paper  develops  a  reconstruction  methodology  for  utilizing  XLT  to  perform 
depth-resolved  imaging  in  a  geometry  appropriate  for  tumor-resection  applications.  This 
method  develops  a  hybrid  X-ray /optical  reconstruction,  which  allows  XLT  spatial 
encoding  in  a  limited  angle  geometry,  and  diffuse  optical  spatial  discrimination  for  the 


remaining  dimensions.  This  technique  augments  that  of  Pratx  et  al,  who  encoded  all 
spatial  dimensions  [9] ;  such  a  technique  is  more  suitable  for  such  applications  as  small- 
animal  imaging.  The  advantage  of  this  new  approach  is  that  enables  XLT  in  surgical 
applications  such  as  breast  or  brain  excision,  and  may  reduce  dose.  The  performance  of 
this  technique  is  examined  in  both  numerical  and  experimental  phantoms  for  various 
object  sizes  and  positions,  within  a  geometry  that  mimics  breast  and  brain  intraoperative 
geometries. 

2,  Methods 

2. 1  Experimental  Setup 

The  equipment  used  for  this  study  consisted  of  an  X-ray  radiation  source  to  excite  the 
phosphors,  and  an  optical  detector  to  sample  the  photon  fluence.  During  acquisition,  the 
radiation  source  is  collimated  into  a  thin  slice  as  described  by  Pratx  et  al.[9]  to  excite  a 
plane  shaped  volume.  An  optical  camera  samples  the  emitted  light.  This  schematic  is 
shown  in  Figure  la.  The  experimental  setup  is  shown  in  Figure  lb. 

The  X-ray  source  used  for  the  measurements  in  this  paper  was  a  50kVp  X-ray 
superficial  unit  (Pantak  Therapax- 1 50,  Elimpex,  AT)  with  a  10cm  exit-diameter  cone 
applicator.  This  cone  was  placed  17.5cm  from  surface  of  the  phantom.  The  beam  was 
collimated  to  1mm  wide  (verified  optically)  by  carefully  positioned  50mm  thick  lead 
bricks.  A  high  sensitivity  EM-CCD  camera  (Pro-EM,  Princeton  Instruments,  NJ)  with  a 
P/1.4  lens  was  positioned  ~20  cm  from  the  surface  of  the  phantom.  This  distance  was 
chosen  to  minimize  X-ray  photon  noise[8].  As  an  alternative,  optically  clear  leaded 
acrylic  or  leaded  glass  could  be  placed  between  the  camera  and  the  sample  and  used  to 
reduce  X-ray  noise  on  the  CCD  and  allow  the  camera  to  be  placed  closer  to  the  sample  to 
collect  more  light. 

2.2  Image  Formation 

2.2.1  X-ray  Luminescence  Emission  Forward  Model 

To  determine  the  concentration  of  nanophosphors,  the  radiation  (X-ray  and  optical)  must 
be  modeled.  The  emission  of  X-ray  excitable  nanophosphors  is  linearly  dependent  on  the 
dose  imparted  to  the  tissue  [8],  d  (units  of  Gy),  the  fractional  efficiency  of  the  phosphor 


in  converting  ionizing  energy  to  optical  emission,  F,  and  the  concentration,  c  (units  of 
mg/ml).  The  luminescent  photon  density  from  the  nanophosphors,  0,  due  to  an  incident 
radiation  beam  is  therefore: 

^  =  rdc  (1) 

Determining  the  total  ionization  energy  imparted  to  the  tissue  (dose)  is  a  procedure  that 
requires  calibration  to  incorporate  the  properties  of  the  radiation  emitted  by  the  X-ray 
system.  This  calibration  is  system-specific,  taking  into  account  tube  potential,  geometry. 
X-ray  tube  target,  and  filter  material;  these  factors  taken  together  form  an  X-ray 
spectrum,  known  as  the  beam  quality.  The  dose  at  depth  is  determined  using 
measurements  from  a  calibrated  ionization  chamber  in  a  phantom  and  composed  into  a 
look-up-table,  the  Percent  Depth  Dose  (PDD)  curves.  This  system-wide  calibration  is 
performed  periodically [20].  Using  the  PDD  curves,  dose  at  a  specific  depth  in  the  tissue 
can  be  determined  by  knowing  the  source-to-surface  distance  between  the  X-ray  tube  and 
the  tissue.  This  method  can  have  high  quantitative  accuracy  of  1-2%  [21].  Another 
method  to  accurately  determine  dose  is  through  Monte-Carlo  methods,  which  model  the 
system,  including  the  above  factors  and  also  including  patient  anatomy.  This  method  can 
calculate  dose  with  high  accuracy  as  long  as  comprehensive  modeling  of  beam  quality  is 
performed[22].  In  this  study,  we  used  the  PDD  curves  to  determine  dose. 

2.2.2  Dijfuse  Optical  Forward  Model 

Images  acquired  at  the  tissue  surface  are  input  into  a  photon  propagation  model  to 
determine  the  phosphor  distribution.  The  images  from  the  CCD  camera  are  first 
processed  to  remove  X-ray  noise  using  a  simple  gradient-threshold  algorithm,  and  then 
input  into  the  algorithm  as  the  data.  The  propagation  of  optical  light  can  be  approximated 
by  the  lossy  photon  Diffusion  Equation  (DE)  [23],  which  yields  the  photon  density  in 
tissue.  The  DE  is  valid  for  many  soft  human  tissues,  including  the  breast,  lung,  prostate, 
brain,  etc.  [24] .  Eollowing  excitation  from  X-ray  radiation,  the  time-independent 
luminescence  photon  density  emitted  from  the  nanophosphors  is: 

0(r)  =  -V-  D(r)V0(r)  + At,(r)<^(r)  (2) 

where  ^r)  is  the  photon  fluence  at  position  r,  in  units  of  photons  per  area  per  time,  and 
0(r)  is  the  photon  density,  in  units  of  photons  per  volume  per  time.  Photon  propagation 


is  affected  by  the  absorption  and  diffusion  coefficients  of  the  tissue,  /r„  and  D, 
respectively,  which  are  dependent  on  wavelength.  The  diffusion  coefficient,  D,  is 

defined  as  D  = - ^ - ,  where  /r ,  ’  is  the  reduced  scattering  coefficient  of  the  emitted 

photons.  A  type  III  boundary  condition  is  used  to  model  the  photon  fluence  at  the 
boundary,  -DS/cj)-  n  =  a(j),  where  a  defines  the  internal  reflection  of  the  light  at  the  tissue 
boundary  due  to  the  index  of  refraction  mismatch  between  tissue  and  air[25][26],  and  the 
unit  vector  n  is  normal  to  the  surface  of  the  phantom.  Because  no  unique  solution  exists 
for  (2)  with  arbitrary  boundaries,  equation  (2)  is  approximated  with  the  Finite  Element 
Method(FEM)[23].  This  problem  is  similar  to  the  diffuse  optical  fluorescence  model 
introduced  by  Jiang[27],  and  is  adapted  here. 

As  described  by  Jiang,  the  photon  emission  may  be  approximated  with  the  FEM 
by: 

[AM  =  {b}  (3) 

where  A  is  the  FEM  approximation  of  the  physics  of  photon  propagation  (the  right-hand- 
side  of  equation  (2))  and  b  is  the  approximation  to  the  light  source  (the  left-hand-side  of 
equation  (2)).  More  specifically,  the  physics  of  the  photon  propagation  is  approximated 
with  the  FEM  by: 

Aij  =  {-DA/ipj  -A/ip.-  (4) 

where  are  the  volume  elements  that  discretize  the  imaging  domain  and  form  a 
geometrical  mesh  defined  over  the  entire  imaging  domain.  A  is  integrated  over  this 
imaging  domain.  The  source  (in  this  case,  the  light  emitted  from  the  phosphors  which 
were  excited  by  the  X-ray  source)  and  boundary  integral  are  approximated  with  the  FEM 
by: 

In  \  m 

\j=l  /  j  =1  boundary 

The  time  component  of  the  luminescence  lifetime  is  ignored  since  the  measurements  in 
this  work  are  from  an  integrating  CCD  camera,  and  the  measurement  time  is  much 
greater  than  the  luminescence  lifetime;  effects  from  the  minimal  afterglow  of  the 
phosphors  are  ignored. 


The  FEM  model,  G,  is  used  to  generate  estimates  for  the  photon  fluence 
given  the  optical  properties  of  the  tissue,  the  concentration  of  phosphors,  c,  and  the  FEM 
mesh.  An  estimate  for  the  photon  fluence,  (p,  can  be  calculated  by  solving  equation  (3): 

(t>  =  G  =  [AY\b}  (6) 

In  this  paper,  the  imaging  domain  is  known,  and  is  assumed  that  the  endogenous  optical 
tissue  properties  are  known,  so  the  model  is  dependent  only  on  the  unknown,  c.  Figure  2 
shows  the  optical  photon  fluence  for  a  numerical  phantom  with  100:1  phosphor 
concentration  between  an  object  and  the  background.  Two  different  source 
configurations  are  shown,  each  with  the  X-ray  direction  of  propagation  in  the  horizontal 
(left/right)  direction.  In  Figure  2a-b,  the  X-ray  source,  indicated  by  the  red  circle, 
irradiates  a  horizontal  line  passing  through  the  background,  whereas  in  Figure  2c-d,  the 
X-ray  irradiates  the  horizontal  line  passing  through  the  middle  of  the  phosphor- 
containing  object. 

2.2.3  X-ray  Nanophosphor  Concentration  Reconstruction 

The  goal  of  X-ray  Fuminescence  Tomography  is  to  determine  the  phosphor  distribution. 
This  spatial  distribution  can  be  determined  by  minimizing  the  difference  between  the 
measured  photon  flux  from  the  camera,  (^^(r),  and  the  simulated  photon  flux,  at 

identical  sample  locations.  This  is  accomplished  by  minimizing  the  F2-norm  of  the 
objective  function  in  an  optimization  routine: 

(7) 

where  Q  is  the  objective  function  to  minimize.  Because  of  the  large  dynamic  range,  the 
first  term  in  (7)  is  formulated  from  the  natural  logs  of  the  photon  fluxes.  This  is  an 
underdetermined  problem,  as  measurements  are  made  only  at  the  boundary.  Because  this 
problem  is  underdetermined,  the  model,  G,  is  linearized  with  a  Taylor  approximation  and 
formed  into  an  iterative  algorithm  as: 

G(c,)  =  G(c,.,)  +  G'Ac 

where  G’  is  the  partial  differential  of  the  model  with  respect  to  the  concentration,  also 
known  as  the  Jacobian,  J.  Minimizing  equation  (7)  with  respect  to  c  and  substituting 
G(cJ  from  equation  (8)  into  (j)^  and  /  for  G’  yields: 


(9) 


2/((G(C;.i)  +  /Ac)  -  ^^ )  =  0 
Solving  for  the  concentration  yields: 

Ac  =  -[/V]V(0,-0^) 

This  problem  is  ill-posed,  so  it  is  solved  using  the  Levenberg-Marquardt  [28]  algorithm, 
which  includes  a  stabilization  parameter.  A,  in  the  inversion  to  avoid  singularities: 

Ac  =  -[/V+A/]V((^,-<^^) 

Equation  (5)  is  iterated  until  a  minima  is  reached  (the  L2  norm  of  the  update  is  less  than 
1%  of  the  previous  iteration),  or  until  15  iterations  are  performed,  whichever  occurs 
earlier.  The  stabilization  parameter  is  reduced  at  each  iteration  as  the  algorithm 
approaches  the  minimum  and  converges  on  the  solution. 

Although  this  study  focused  on  applications  where  a  single  angle  is  ideal,  note 
that  this  algorithm  is  not  limited  to  a  single  angle.  Thus,  this  algorithm  is  appropriate  for 
any  sparse-angle  geometry. 

2.3  Phantom  Study 

The  performance  of  the  experimental  setup  and  the  reconstruction  algorithm  were  tested 
by  varying  the  size  and  location  of  a  lesion-simulating  object.  The  relationships  between 
source-object  and  detector-object  distance  on  resolving  an  object  of  various  sizes  were 
determined  with  both  numerical  and  experimental  phantoms.  The  metrics  used  to 
determine  system  performance  were  object  location  and  object  size.  Location  error  in 
both  the  X-ray  excitation  and  optical  read-out  dimensions  was  determined  by  calculating 
the  distance  between  the  true  centroid  of  the  lesion  and  the  location  of  the  maximum 
recovered  value  of  the  phantom.  Object  size  in  both  the  X-ray  excitation  and  optical 
read-out  dimensions  were  determined  by  calculating  the  full  width  at  half  maximum 
(FWHM)  of  the  object  in  these  dimensions.  Concentration  sensitivity  and  contrast 
recovery  were  also  examined. 

2.3.1  Numerical  Phantoms 

Numerical  phantoms  were  utilized  to  test  the  position-accuracy,  object  size-accuracy,  and 
sensitivity  of  the  algorithms.  The  position  and  object-size  phantoms  investigated  the 


recovery  of  simulated  tumors  with  10: 1  contrast  between  the  object  and  the  background 
and  phosphor  concentrations  of  lOpg/ml.  These  lesions  varied  in  size  between  2-14mm, 
and  were  placed  at  different  locations  in  the  phantom,  as  depicted  in  Figures  3a  and  4a 
(note  the  multiple  arrows,  which  indicate  the  minimum  and  maximum  extent  of  the  object 
locations  investigated).  Sufficient  dose  (IcGy)  was  given  to  yield  signal-to-noise  (SNR) 
greater  than  10  for  each  object  position  -  this  methodology  allowed  a  performance  test  of 
the  algorithm  for  all  object  positions. 

The  sensitivity  phantom  included  an  object  with  varying  concentration  (Figure 
5a),  and  varying  contrast  (Figure  5b).  The  6mm-diameter  object  was  placed  at  the  center 
of  the  phantom  along  the  dimension  of  the  detectors  (the  long-axis),  and  moved  at 
various  depths  away  from  the  detectors.  We  used  the  phosphor  properties  from 
Kandarakis  et  al.  [29]  to  obtain  quantification  of  the  emitted  light  efficiency  for  their 
lanthanum  oxysulfide:terbium  phosphor,  which  was  1.39x10^^  optical  photons  /  (Gy  x 
mg).  We  incorporated  solid-angle  losses  as  well  as  losses  due  to  lens  inefficiency  (DO- 
1795,  Navitar  Imaging  Solutions,  Rochester,  NY).  SNR  below  10  was  assumed  to  be  too 
low  to  detect. 

All  phantoms  were  2-dimensional,  and  measured  12cm  x  6cm.  Detectors  were 
placed  along  the  long-axis  of  the  phantom,  while  the  collimated  X-ray  source,  1mm  wide, 
was  scanned  along  the  short  axis.  The  phosphor  used  for  this  experiment  mimicked 
GOS:Eu,  demonstrated  in  X-ray  luminescence  imaging  in  a  previous  study,  which  has  a 
strong  luminescence  emission  at  ~618nm.  Background  optical  properties  were  similar  to 
that  of  breast  tissue  [30]  (p^  =  0.0027,  p/  =  0.717). 

2.3.2  Experimental  Phantom 

The  experimental  phantom  is  shown  in  Figure  6.  Figure  6a-b  show  the  relative  layout  of 
the  objects  in  the  phantom,  while  Figure  6c  shows  the  camera-eye  view  of  the  phantom. 
Figure  6d  shows  an  image  of  the  phantom  while  the  phantom  is  irradiated  by  the  X-ray 
beam.  The  optically  clear  acrylic  phantom,  measuring  12cm  x  6cm,  was  filled  with  india 
ink  to  mimic  optical  absorption,  and  intralipid  to  mimic  optical  scatter,  at  the  appropriate 
concentrations.  The  optical  properties  were:  (p^  =  0.0027,  p,’  =  0.717,  as  determined 
from  a  diffuse  optical  spectroscopy  system).  GOS:Eu  phosphor  at  a  concentration  of 


lOmg/ml  was  added  to  two  cylindrical  inclusions,  one  5mm  in  diameter,  and  one  10mm 
in  diameter,  which  were  located  3cm  and  9cm  from  the  edge  nearest  the  X-ray  source, 
respectively.  The  inclusions  were  both  imaged  at  various  depths  from  the  edge  nearest 
the  detector:  10mm,  15mm,  20mm,  30mm;  these  dimensions  are  shown  more  clearly  in 
Figure  7a.  The  exposure  times  and  gains  were  1.5s  at  gain  500,  3.5s  at  gain  800,  3s  at 
gain  1000,  and  7.5s  at  gain  1000,  for  increasing  depth.  Dose  to  the  phantom  varied 
depending  on  the  phosphor  depth,  so  that  a  high  SNR  could  be  acquired  while  the 
phosphors  were  irradiated.  The  doses  to  the  phantom  were  6.7  cGy,  15.6  cGy,  13.4  cGy, 
and  33.4  cGy,  for  increasing  depth. 

3,  Results 

3. 1  Numerical  phantom  results 

As  described  above,  a  lesion-simulating  object  was  placed  in  various  locations  in  the 
volume  so  that  the  performance  of  LAXLT  could  be  analyzed.  The  lesion  location  was 
varied  with  respect  to  the  source-axis  and  detection-axis  separately  to  determine  the 
effects  of  source-object  distance  and  detector-object  distance  on  the  ability  to  resolve  the 
object.  Figure  3  shows  the  results  of  maintaining  a  fixed  detector-object  distance  while 
varying  the  source-object  distance  (varying  depth  with  respect  to  the  X-ray).  These 
dimensions  are  depicted  in  Figure  3a.  The  location  error  for  the  2- 14mm  objects  is 
shown  in  Figure  3b.  This  result  demonstrates  the  excellent  ability  of  the  algorithm  to 
spatially  resolve  an  object,  as  all  errors  in  distance  are  lower  than  3.5mm.  Figure  3c 
shows  the  recovered  object  size  in  the  direction  of  the  X-ray  excitation,  which 
demonstrates  the  insensitivity  with  respect  to  depth,  with  an  average  standard  deviation  in 
the  error  in  recovered  size  of  0.82mm  for  varying  source-object  distances.  In  Figure  3d, 
the  recovered  size  of  the  object  in  the  optical  dimension  is  examined  as  a  function  of 
source-object  depth.  The  algorithm  is  able  to  distinguish  the  varying  sizes  of  the  objects, 
with  a  standard  deviation  in  the  error  in  recovered  size  of  0.83mm  for  the  varying  source- 
object  distance.  There  is  a  slight  tendency  for  blurring  in  the  optical  dimension  due  to  the 
scattering  of  the  optical  photons. 

Figure  4  shows  the  effects  of  varying  the  depth  of  an  object  with  respect  to  the 
optical  dimension.  Figure  4a  shows  the  physical  dimensions  of  the  phantom.  The  object 


was  moved  from  0.5cm  to  4.5cm  from  the  detector,  while  remaining  at  a  fixed  distance 
from  the  X-ray  source.  Figure  4b  shows  the  recovered  location  error  of  the  centroid  of 
the  object  with  respect  to  depth  from  the  optical  detector.  In  this  experiment,  the  ability 
to  resolve  an  object  of  ~2mm  was  at  the  limit  of  the  system  with  the  1mm  collimated  X- 
ray  beam  that  was  used.  However,  the  objects  sized  4mm  and  larger  were  resolved  with 
higher  accuracy.  As  expected,  location  error  increases  slightly  as  the  depth  increases,  and 
larger  objects  are  resolved  more  accurately  than  smaller  objects.  The  advantage  of  using 
the  X-ray  is  highlighted  here,  as  even  at  4.5cm  in  depth,  the  object  is  resolved  with  less 
than  1cm  total  error.  Figure  4c  shows  the  recovery  of  object  size  in  the  X-ray  dimension, 
similar  to  Figure  3c.  The  mean  standard  deviation  for  the  error  in  recovered  object  size 
was  1.36mm.  Figure  4d  shows  that  the  ability  to  determine  the  size  of  the  object  in  the 
optical  dimension  has  increased  variability  compared  to  the  X-ray  dimension,  as  the 
mean  standard  deviation  of  the  error  in  recovered  size  with  respect  to  depth  is  2.4mm. 
This  larger  standard  deviation  is  consistent  with  object-blurring  at  depth  from  the  optical 
detector;  this  is  expected  due  to  the  ill-posed  nature  of  the  algorithms  needed  for  the 
optical  photon  modeling.  Similar  to  Figure  3d,  this  system  is  able  to  distinguish  between 
the  different  sized  objects  in  the  optical  dimension,  although  again,  there  is  a  tendency  for 
dilation  in  the  optically-resolved  dimension. 

Overall,  the  numerical  phantom  results  show  localization  accuracy  with  median 
error  of  2.2mm  (mean  of  6.3mm),  or  4.1%  (mean  of  1 1.5%)  of  object  depth  for  all 
lesions.  The  recovered  lesion  size  has  lower  size-bias,  with  median  error  of  -8.1%  vs. 
87.5%  (mean  of  19.4%  and  1 18.3%,  respectively),  in  the  X-ray  excitation  direction  vs. 
the  optical  direction,  respectively.  Again,  this  optical  dilation  is  expected  due  to  the 
increased  optical  scatter  compared  to  X-ray.  This  technique  is  invariant  in  recovered  size 
with  respect  to  depth,  as  the  standard  deviation  is  less  than  2.5mm. 

The  concentration  phantom  is  shown  in  Figure  5  a.  The  required  dose  (in  Gy)  to 
reach  an  SNR  of  10  is  plotted  for  varying  concentrations.  It  is  apparent  from  this 
calculation  that  a  6mm-diameter  object  with  pg/ml  concentration  is  detectable  in  this 
geometry  with  standard  CT  doses.  With  the  doses  currently  used  in  lORT, 
concentrations  to  ng/ml  are  detectable  at  depth.  As  shown  in  Figure  5b,  with  sufficient 
signal,  contrast  can  be  recovered  for  all  contrast  to  background  ratios  tested,  for  depths  up 


to  45mm.  Here,  the  advantage  of  the  ability  of  the  collimated  X-ray  to  selectively  excite 
the  phosphors  is  clear. 

3.2  Experimental  phantom  results 

Photographs  of  the  experimental  phantom  are  shown  in  Figure  6;  the  schematic  of  the 
experimental  phantom  is  shown  in  Figure  7a,  and  the  results  are  shown  in  Figure  7b-d. 
Reconstructions  for  each  case  are  shown  in  Figure  8,  with  increasing  depth  with  each 
row.  In  this  experiment,  a  1-  5mm  object,  and  a  1-  10mm  object  were  imaged  at  various 
depths.  For  all  depths,  the  recovered  location  error  for  both  objects  is  less  than  6mm,  and 
is  independent  of  object  size.  Positional  errors  are  recovered  within  8.6%  of  the  effective 
depth  for  a  5mm  object,  and  within  5.2%  of  the  depth  for  a  10mm  object.  Similar  to  the 
results  from  the  numerical  phantom,  the  ability  to  discern  object  size  is  highly  accurate  in 
the  X-ray  dimension,  and  is  invariant  with  depth.  The  ability  to  resolve  the  object  in  the 
optical  dimension  is  accurate  to  within  2mm.  Object  size  median  error  is  within  2.3% 
and  2%  for  the  5mm  and  10mm  objects,  respectively. 

4,  Discussion 

Molecular  imaging  has  long  been  identified  as  potentially  a  vital  tool  in  surgery  to  aid  in 
the  identification  of  important  tissue  structures  and  tumor  tissue[31].  Until  recently,  the 
use  of  molecular  imaging  in  surgery  has  been  limited,  both  due  to  the  lack  of  appropriate 
imaging  tools,  as  well  as  the  lack  of  highly-specific  contrast  agents  (the  problem  with 
most  endogenous  fluorescence  agents  in  the  body).  Recent  developments  in  both  of  these 
areas  have  increased  the  interest  of  molecular  surgical-guidance  [3].  For  instance. 
Fluorescence  Molecular  Tomography  (FMT)  shows  great  promise  in  positively  affecting 
surgical  outcomes,  as  it  is  non-ionizing,  has  a  wealth  of  knowledge  of  contrast  agents, 
and  can  be  implemented  at  low  cost.  FMT  has  found  utility  in  tumor-margin  excision 
[32],  sentinel  lymph-node  mapping  [4],  and  avoidance  of  critical  structures  such  as 
nerves[33].  Alternatively,  sentinel  lymph  node  identification  with  radio-emission 
imaging  is  widely  adopted[5];  its  use  in  surgical  resection  is  not  well  established. 

The  lack  of  wide-adoption  of  molecular  imaging  in  surgical  applications  can 
partly  be  blamed  on  the  limited  flexibility  of  these  imaging  technologies.  While  FMT  is 


ideal  for  structures  at  the  surface  and  several  millimeters  beyond,  its  limitations  in  depth 
penetration  preclude  its  ability  to  unmask  disease  that  exists  several  centimeters  beneath 
the  surface.  Radio-emission  imaging  is  limited  because  of  the  need  to  image  gamma 
photons  which  are  highly  penetrating  in  tissue  -  this  prevents  the  identification  of  lesion 
depth,  and  increases  the  burden  on  the  surgeon  to  patiently  dig  through  the  tissue  until  the 
object  is  identified.  In  the  case  of  a  sub-millimeter  sized  tumor  tissue,  this  practice  would 
be  overly  laborious. 

Clearly,  a  void  exits  in  the  ability  of  molecular  instrumentation  to  resolve 
millimeter  or  sub-millimeter  objects  at  depths  of  several  centimeters.  LAXLT  may  be 
able  to  fill  this  void.  Its  advantage  is  the  ability  to  resolve  objects  at  several  centimeters 
of  deep.  Figure  9  demonstrates  this  ability,  as  it  is  shown  that  both  the  5mm  and  10mm 
objects  in  the  experimental  phantom  were  able  to  be  resolved  at  a  detector-object  depth  of 
beyond  3cm.  This  is  a  stark  contrast  to  FMT,  as  Kepshire  et  al.  [18]  demonstrated  that 
depth  linearity  degrades  beyond  1cm.  If  the  depth  penetration  of  the  X-ray  source  is 
considered,  LAXLT  is  linear  beyond  10cm.  This  increased  performance  is  due  to  the  low 
scatter  and  high  penetrability  of  the  X-ray  excitation,  which  allows  more  ideal  imaging 
geometries  to  be  chosen  for  surgical  guidance. 

The  benefit  of  the  technique  developed  in  this  study  compared  to  previous 
developments  with  XLT  is  three-fold:  increased  temporal  resolution  is  possible  because 
of  the  optical  read-out  of  two-dimensions,  which  eliminates  X-ray  encoding  in  those 
dimensions;  lower  dose  is  possible  because  of  the  decreased  X-ray  excitation;  most 
importantly,  this  technique  is  suitable  for  geometries  where  full  angular  X-ray  encoding 
is  not  possible.  Therefore,  this  technique  enables  XLT  to  be  used  for  surgeries  such  as 
breast  lumpectomy;  we  foresee  applications  such  as  this  to  be  one  of  the  most  important 
future  applications  in  imaging.  In  breast  lumpectomy,  a  full-angular  encoding  with  X-ray 
is  not  desirable  because  deep  critical  structures  are  irradiated.  By  instead  implementing  a 
limited-angle  technique,  only  the  breast  may  be  irradiated.  For  lumpectomy,  this 
technique  would  be  desirable  to  verify  the  position  of  the  lesions  in  the  surgical  supine 
position,  and  to  visualize  remaining  disease  after  resection. 

In  comparing  this  technique  to  other  image-guided  surgical  modalities  such  as 
FMT,  the  advantages  of  LAXLT  are  clear.  LAXLT  was  able  to  resolve  4mm  objects  at 


detector-object  depths  greater  than  4.5cm,  compared  to  FMT  which  is  limited  in  to  about 
1cm  in  resolving  object  dimensions [18].  The  most  impressive  aspect  indicated  by  these 
results  is  the  reduction  in  the  blurring  of  the  object  with  respect  to  depth;  this  highlights 
the  advantage  of  incorporating  the  X-ray  excitation,  which  due  to  its  relatively 
insignificant  scatter  at  these  depths,  can  pinpoint  the  depth  of  the  object  with  respect  to 
the  optical  depth  dimension.  The  optical  read-out  is  then  used  to  determine  the  other 
dimensions.  Since  the  depth  is  known  to  high  accuracy  (shown  in  Figure  4),  the  diffuse 
algorithm  properly  models  the  diffuse  nature  of  the  light,  and  significantly  reduces  object 
blurring,  as  presented  in  Figure  5.  Thus,  with  LAXLT,  dose  and  imaging  time  are 
reduced  significantly  compared  to  XLT. 

It  is  intuitive  that  the  resolution  will  be  limited  by  the  width  of  the  X-ray  beam,  at 
shallower  depths.  At  deeper  depths,  the  scatter  of  the  X-ray  beam  should  be  taken  into 
account.  In  this  study,  a  I  mm  collimated  X-ray  beam  was  able  to  successfully  visualize  a 
2mm  object  up  to  2cm,  and  a  4mm  object  up  to  4.5cm.  Higher  resolution  should  be 
attained  with  a  narrower  beam;  we  are  currently  investigating  this  effect.  The  ability  to 
resolve  an  object  at  depth  is  instead  limited  by  dose. 

Although  this  study  was  adapted  to  geometries  best  suitable  to  tumor  resection,  it 
should  be  noted  that  the  algorithm  presented  in  this  work  is  generalizable  to  any 
geometry;  it  is  especially  useful  for  sparse-angle  geometries.  An  additional  use  for  this 
technique  would  be  for  intraoperative  probes,  where  simple  coregistration  (assuming  the 
catheter  is  radio-opaque)  between  the  X-ray  source  and  the  optical  detection  catheter 
could  provide  assessment  of  molecular  status  at  a  remote  location. 

An  increase  in  dose  to  potentially  healthy  tissue  is  one  disadvantage  of  this 
technique.  Although  we  calculated  that  the  maximum  dose  to  the  phantom  was  33.4  cGy, 
this  dose  will  be  reduced  to  6cGy  with  a  more  favorable  optical  setup  where  the  camera  is 
closer  to  the  object.  This  dose  may  be  reduced  further  with  a  contact  optical  setup.  Still, 
this  dose  may  preclude  its  use  in  a  screening  setting,  especially  due  to  the  increasing 
awareness  of  increasing  radiation  exposure  in  medicine[34].  However,  post-surgical 
radiation  therapy  is  commonly  prescribed  as  a  means  to  destroy  cancer  cells  that  may  not 
have  been  removed  around  the  margin  [35,  36].  In  Accelerated  Partial  Breast  Irradiation, 
5-20  Gy  of  radiation  is  given  to  the  resected  cavity  in  1-5  fractions,  to  reduce  the 


morbidity  of  whole  breast  radiation  therapy  [37].  In  this  context,  LAXLT  may  have  great 
utility  to  identify  larger  regions  that  may  have  been  missed  during  surgery,  and  may 
require  subsequent  surgical  investigation.  Its  sensitivity  at  this  dose  might  enable  the 
identification  of  micro-disease,  or  a  very  low  concentration  to  be  injected. 

The  other  disadvantage  of  this  technique  is  its  use  of  nanoparticles.  Although  we 
have  demonstrated  low-toxicity  of  our  nanoparticles  in  cells,  nanoparticles  will  have 
different  effects  in  a  living  system.  This  topic  is  beyond  the  scope  of  this  paper,  but  it  is 
important  to  recognize  that  cancer  nanotechnology  is  a  major  venture  in  the  National 
Cancer  Institute,  and  this  technique  will  benefit  from  the  knowledge  gained  from  this 
program. 

5,  Conclusions 

Limited  Angle  X-ray  Luminescence  Tomography  has  been  developed  as  a  means  to 
image  molecular  deeper  than  is  available  with  Fluorescence  Molecular  Tomography.  A 
reconstruction  algorithm,  based  on  a  hybrid  X-ray  excitation/diffuse  optical  emission 
model,  was  tested  in  a  numerical  and  experimental  phantom  that  had  dimensions  similar 
to  the  human  breast.  It  was  found  that  objects  as  small  as  2mm  in  diameter  could  be 
resolved  at  depths  of  up  to  4-5cm.  It  was  shown  that  the  depth  of  the  object  with  respect 
to  the  X-ray  source  position  had  little  effect  on  object  recovery  in  this  volume.  It  was 
then  demonstrated  experimentally  that  both  a  5mm  and  10mm  object  could  be  resolved  at 
depths  of  at  least  3cm  -  these  results  agreed  with  the  numerical  phantom  results,  thus 
validating  the  simulation.  If  the  challenges  to  engineering  biocompatible  phosphors  can 
be  resolved,  LAXLT  may  have  utility  in  surgical  applications  where  small  lesions  must 
be  imaged  at  depth  of  a  few  centimeters,  such  as  during  breast  lumpectomy  surgery. 
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Figure  Captions 


Figure  1 :  (a)  Schematic  of  LAXLT  while  imaging  a  phantom  with  a  single  spherical 
object  (b)  Experimental  setup  used  in  this  study.  A  black-box  was  used  to  enclose  the 
phantom  (not  shown)  and  eliminate  ambient  light  from  the  experiment. 

Figure  2:  (a)  Numerical  phantom  with  a  single  object  located  at  the  center  of  the  domain, 
and  a  10: 1  contrast  in  phosphor  between  the  object  and  background.  The  X-ray  source 
irradiates  the  domain  along  the  long-axis  of  the  phantom  at  a  position  of  X=10,  Y=- 
60:60.  (b)  The  emission  fluence  field  of  (a),  (c)  Numerical  phantom  with  the  X-ray 
source  irradiating  at  a  position  of  X=40,  Y=-60:60.  (d)  The  emission  fluence  field  of  (c). 

Figure  3:  (a)  Numerical  phantom  experiment  examining  the  effect  of  moving  a  variable 
sized  object  away  (4  positions,  indicated  by  the  arrows)  from  the  X-ray  source  plane  (S), 
while  keeping  the  distance  from  the  detection  plane  (D)  fixed,  (b)  Recovered  object 
location  error,  (c)  Recovered  object  diameter  FWHM  with  respect  to  the  X-ray  source 
dimension,  (d)  Recovered  object  diameter  FWHM  with  respect  to  the  optical  detection 
dimension. 

Figure  4:  (a)  Numerical  phantom  experiment  examining  the  effect  of  moving  a  variable 
sized  object  away  (8  positions,  indicated  by  the  arrows)  from  the  optical  detection  plane 
(D),  while  keeping  the  distance  from  the  X-ray  source  plane  (S)  fixed,  (b)  Recovered 
object  location  error,  (c)  Recovered  object  diameter  FWHM  with  respect  to  the  X-ray 
source  dimension,  (d)  Recovered  object  diameter  FWHM  with  respect  to  the  optical 
detection  dimension. 

Figure  5:  (a)  Numerical  phantom  experiment  examining  the  effect  of:  (a)  variable  object 
to  background  contrast  and  (b)  variable  concentration,  of  a  6mm-diameter  object  vs. 
depth  from  the  detection  plane. 


Figure  6:  Experimental  phantom:  (a-b)  Two  views  from  above  the  black-box  including 
the  phantom  and  collimator  bricks,  (c)  Image  taken  by  the  high-sensitivity  CCD  camera 
of  the  phantom  under  ambient  light  and  no  irradiation,  (d)  CCD  image  with  no  ambient 
light  and  X-ray  irradiation. 

Figure  7:  (a)  Experimental  phantom  of  two  objects  moved  away  from  the  optical 
detection  plane  (D)  at  4  positions,  (b)  Recovered  object  location  error,  (c)  Recovered 
object  diameter  FWHM  with  respect  to  the  X-ray  source  dimension,  (d)  Recovered  object 
diameter  FWHM  with  respect  to  the  optical  detection  dimension. 

Figure  8:  Image  reconstructions  of  2  objects  with  depths  of  (a)  lOmm,  (b)  I5mm,  (c) 
20mm,  and  (d)  30mm  from  the  detection  plane. 


Figure  9:  (a)  Object  location  linearity  of  the  depth  of  the  object  with  respect  to  the 
detection  plane,  (b)  Object  location  linearity  of  both  the  depth  of  the  object  with  respect 
to  the  detection  plane  summed  with  distance  with  respect  to  the  source  plane. 
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Radioluminescent  nanophosphors  enable 
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Abstraet:  We  demonstrate  the  ability  to  image  multiple  nanoparticle-based 
contrast  agents  simultaneously  using  a  nanophosphor  platform  excited  by 
either  radiopharmaceutical  or  X-ray  irradiation.  These  radioluminescent 
nanoparticles  emit  optical  light  at  unique  wavelengths  depending  on  their 
lanthanide  dopant,  enabling  multiplexed  imaging.  This  study  demonstrates 
the  separation  of  two  distinct  nanophosphor  contrast  agents  in  gelatin 
phantoms  with  a  recovered  phosphor  separation  correlation  of  -0.98.  The 
ability  to  distinguish  the  two  nanophosphors  and  a  Cerenkov  component  is 
then  demonstrated  in  a  small  animal  phantom.  Combined  with  the  high- 
resolution  potential  of  low-scattering  X-ray  excitation,  this  imaging 
technique  may  be  a  promising  method  to  probe  molecular  processes  in 
living  organisms. 

©2012  Optical  Society  of  America 

OCIS  codes:  (1 10.4234)  Multispectral  and  hyperspectral  imaging;  (170.0170)  Medical  optics 
and  biotechnology;  (160.4236)  Nanomaterials. 
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1.  Introduction 

The  application  of  luminescent  nanoparticles,  such  as  quantum  dots,  in  diagnostic  molecular 
imaging  is  receiving  much  attention  due  to  the  potential  to  improve  the  differentiation 
between  normal  and  malignant  tissue  by  identifying  signatures  of  disease  [1].  Eor  example, 
multicolored  quantum  dots  enable  the  ability  to  probe  multiple  signatures  of  disease 
simultaneously;  this  is  recognized  to  be  potentially  more  beneficial  than  visualizing  a  single 
species  [2].  Examples  of  these  molecular  markers  include  the  Uvps  integrin  and  VEGF 
receptors  as  indicators  for  vascular  angiogenesis  [3],  and  the  translocation  of 
phosphatidylserine  for  apoptosis  [4].  The  advantage  of  nanoparticles  in  this  context  is  that 
they  may  be  tuned  to  emit  at  varying  wavelengths,  without  severely  altering  their  physical 
structure  [5].  Although  quantum  dots  have  been  investigated  extensively  for  biological 
applications  [6,  7],  nanophosphors  have  been  relatively  ignored.  Like  quantum  dots, 
nanophosphors  have  high  quantum  efficiencies,  sharp  emission  peaks,  high  photostability,  and 
emission  wavelength  tunability  [8,  9].  Nanophosphors  have  a  distinct  advantage  in  that  they 
may  emit  light  when  irradiated  by  ionizing  radiation,  a  process  called  radioluminescence. 
This  work  demonstrates  the  multiplexing  potential  of  these  nanophosphors,  using  an 
exemplary  nanophosphor  system. 

While  radioluminescence  of  phosphor  material  has  long  been  used  in  radiation  detectors, 
the  use  of  radioluminescent  nanophosphors  (RLNPs)  in  biological  contexts  is  just  beginning 
to  be  explored  [10].  Novel  applications  currently  under  investigation  for  radioluminescent 
imaging  (RLI)  include  direct  biological  feedback  of  molecular  changes  and  tumor  burden 
during  radiation  treatment.  This  information  may  be  used  as  a  means  to  enhance  radiation 
treatment  efficacy.  In  addition,  RLI  may  be  used  to  probe  molecular  information  at  a  superior 
resolution  to  other  molecular  imaging  modalities,  such  as  Fluorescence  Molecular 
Tomography  (FMT),  Positron  Emission  Tomography  (PET),  and  Single  Photon  Computed 
Emission  Tomography  (SPECT),  by  collimating  the  excitation  X-ray  beam.  This  high- 
resolution  technique,  called  X-ray  Luminescence  (Computed)  Tomography  (XLCT),  utilizes 
narrowly-collimated  X-ray  beams  to  resolve  deep-seated  RLNPs.  This  is  made  possible 
because  of  the  very  low  X-ray  scatter  in  small-animals  [11,  12]. 

This  study  demonstrates  the  potential  use  of  RLNPs  doped  with  different  rare-earth 
luminescent  centers  to  enable  multiplexed  RLI  after  excitation  by  ionizing  radiation.  In 


particular,  it  demonstrates  the  ability  to  recover  the  concentrations  of  each  RLNP  from 
multispectral  images  of  the  subject.  To  highlight  the  flexibility  of  this  technique,  multiplexed 
imaging  was  demonstrated  using  two  different  excitation  schemes:  radiopharmaceuticals  (**F) 
and  X-ray  radiation. 

2.  Methods 

2. 1  Imaging  hardware 

A  custom  imaging  system  was  built  to  image  and  measure  the  concentrations  of  the 
nanophosphor  contrast  agents.  A  512x512  pixel  back-illuminated  deep-cooled  (-70C)  CCD 
camera  (ProEM,  Princeton  Instruments,  Trenton,  NJ)  was  outfitted  with  a  F/#0.95  imaging 
lens  (DO-5095,  Navitar,  Rochester,  NY),  and  leaded  glass  was  placed  above  the  subject. 
Optical  bandpass  filters  (546nm  +/-  lOnm  for  Tb^"^  doped  RLNPs,  700nm  +/-  20nm  for  Eu’"^ 
doped  RLNPs;  Andover  Corporation,  Salem,  NH)  were  placed  in  front  of  the  camera  lens  to 
selectively  collect  emitted  light  in  the  unique  spectral  regions.  A  motorized  platform  adjusted 
the  working  distance  of  the  camera  to  the  subject  and  controlled  the  field  of  view.  This 
imaging  box  was  placed  in  the  path  of  a  superficial  X-ray  unit  (Pantax  Therapax,  East  Flaven, 
CT),  with  voltage  and  current  set  to  the  X-ray  tube  90keV  and  10mA.  This  system  was 
operated  remotely  from  the  control  room  of  the  X-ray  unit,  outside  the  X-ray  room  [11]. 

Optical  spectra  from  the  RLNPs  (Fig.  1(a)  &  Fig.  1  (b))  was  collected  by  measuring  the 
emission  of  dry  powder  samples  in  cuvettes.  The  distal  end  of  a  10m  400pm  optical  fiber  was 
placed  in  contact  with  the  side  of  cuvettes  containing  each  sample,  while  the  proximal  end 
was  attached  to  a  spectrophotometer  (SP2150,  Princeton  Instruments)  connected  to  the  CCD 
camera,  which  was  operated  from  the  X-ray  control  room.  Optical  emission  spectra  within 
the  visible  and  NIR  range  was  acquired  with  custom  software  written  in  Labview  (National 
Instruments,  Austin,  TX).  By  acquiring  vertical  spectral  bins  for  each  wavelength,  image 
processing  was  used  to  remove  X-ray  noise,  which  stochastically  appeared  in  the  CCD. 
Alternatively,  Cerenkov  spectra  was  recorded  by  fitting  a  quadratic  curve  to  data  collected 
from  a  commercial  imaging  system  (IVIS  Spectrum,  Perkin  Elmer,  Waltham,  MA).  The 
signal  emitted  from  the  radiopharmaceuticals  was  verified  via  a  small  animal  Positron 
Emission  Tomography  (PET)  system  (Siemens  R4  Micro-PET,  Malvern,  PA). 
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Fig  1.  (a)  Radioluminecent  spectrum  of  the  Bao.sjYo.jFiiTb^*  particles,  which  emit  primarily  in  the  green,  (b) 
Radioluminecent  spectrum  of  the  Bao.sjYo.sF^iEu^*  particles,  which  emit  primarily  in  the  red  and  near-infrared,  (c) 
Transmission  electron  microscopy  image  of  the  RLNPs.  (d)  Flydrodynamic  size  distribution  of  the  RLNPs  as 
determined  by  dynamic  light  scattering. 


2.2  Contrast  agent  synthesis 


In  this  study,  we  demonstrate  two  possible  dopants,  europium  (Eu^"^)  and  terbium  (Tb^^).  A 
detailed  description  of  the  synthesis  of  these  particles  was  recently  described  by  Sun  et  al. 
[13].  The  host  barium  yttrium  fluoride  (Bao,55Yo3F2)  nanocrystals  were  doped  with  either 
terbium  (0.5%)  or  europium  (0.5%)  to  produce  emission  maxima  at  541nm  and  586nm, 
respectively.  Both  RLNPs  emitted  additional  spectral  peaks  towards  the  near  infrared  (NIR) 
at  629nm  and  692nm,  respectively,  as  seen  from  Fig.  1(a)  &.  Fig.  1(b);  these  emission 
wavelengths  enable  imaging  of  deeper  structures.  The  unique  signatures  emitted  by  each 
phosphor  enabled  the  ability  to  perform  multiplexing.  The  surfaces  of  the  RLNPs  were 
modified  with  poly(ethylene  glycol)  (PEG,  MW  600)  through  a  ligand  exchange  process  to 
achieve  colloidal  stability  in  aqueous  solutions.  These  phosphors  can  be  functionalized  with  a 
biotin  linker  bound  to  the  PEG.  Shown  in  Fig.  1(c)  is  the  transmission  electron  microscopy 
(TEM)  image  of  this  RLNP  platform  with  a  cubic  morphology  and  mean  particle  size  of 
14nm.  The  hydrodynamic  size  was  determined  by  dynamic  light  scattering  (DLS);  particles 
with  hydrodynamic  size  of  ~27nm  were  the  most  common,  as  shown  in  Fig.  1(d). 

2.3  Image  multiplexing 

To  form  images  of  the  concentrations,  c,  of  each  of  these  particles,  the  light  emitted,  (j),  was 
recorded  over  the  several  optical  spectral  regions  defined  by  filter  fi  and  filter  f2  to  form 
images  of  the  subject,  ij  =(p(fi)  and  i2  =<l)(fi).  These  images,  and  the  pre-recorded  reference 
spectra,  e,  for  Bao.55Yo.3F2:Tb^'^,  (er*),  Bao.55Yo.3F2:Eu^'^,  {Eeu),  and  the  Cerenkov  luminescence, 
(Ech),  were  input  into  a  linear-least  squares  algorithm  to  extract  the  contributions  from  each 
nanophosphor.  Reference  spectra  was  collected  with  a  calibrated  spectrometer,  thus  enabling 
a  quantitative  comparison  between  the  nanophosphors.  Since  Cerenkov  light  was  generated 
from  the  F-18  radiotracer,  it  was  necessary  to  account  for  the  Cerenkov  component  for  the 
experiments  involving  radiotracer  excitation.  For  experiments  involving  X-ray  excitation,  the 
Cerenkov  component  was  not  present,  since  Cerenkov  light  arises  from  photons  emitted  from 
highly-energetic  charged  particles  that  are  greater  than  the  energy  of  charged  particles  created 
during  X-ray  excitation. 

2.4  Experimental  demonstration 

To  validate  the  ability  of  this  method  to  image  multiple  nanophosphors  simultaneously,  five 
phantoms  (99%  DI  water,  1%  agarose)  were  fabricated  containing  Tb^"^  and  Eu^"^  doped 
RLNPs  of  linearly  varying  concentrations  (in  units  of  mg/ml)  and  mixed  in  the  following 
ratios:  (0:10;  2. 5:7. 5;  5:5;  7. 5:2. 5;  and  10:0).  To  demonstrate  this  method  in  a  pre-clinical 
mouse  model,  four  batches  of  RLNPs  were  each  mixed  with  52  pCi  (Curie)  of  activity  of  18-F 
and  Matrigel  (BD  Biosciences,  Sparks,  MD)  and  injected  into  both  forelegs  and  flanks  of  a 
nude  mouse,  with  Eu^'^-doped  RLNPs  (Eu)  only  on  the  left  foreleg,  Tb^^-doped  RLNPs  (Tb) 
only  on  the  left  flank,  an  equal  mixture  of  Eu^'^-doped  RLNPs  and  Tb^^-doped  RLNPs  (Mix) 
on  the  right  flank,  and  an  inactive  undoped  nanophosphor  control  (Ctrl)  on  the  right  foreleg. 
For  the  radioisotope-excited  experiment,  the  subject  was  imaged  immediately  after  injection. 
For  the  X-ray-excited  experiment  several  days  later  (after  the  18-F  had  fully  decayed),  the 
same  subject  was  placed  in  the  imaging  apparatus  and  excited  by  a  90keV  X-ray  source, 
which  was  transmitted  through  the  dark-box.  The  same  subject  was  imaged  with  the  X-ray 
and  radiopharmaceutical  excitation  on  separate  days  to  allow  the  radioisotope  to  decay 


thoroughly.  Concentration  images  of  the  animals  were  normalized  due  to  the  lack  of  recorded 
tissue  optical  properties. 

3.  Results 

3.1  Phantom  validation 


Fig.  2(a)  displays  the  unmixed  images  for  linearly  increasing  (top-to-bottom)  amounts  of 
Tb^'^-doped  RLNPs  on  the  left,  and  decreasing  (top-to-bottom)  Eu^^-doped  RLNPs  on  the 
right,  at  each  concentration.  The  median  values  for  each  region  of  interest  were  plotted  with 
respect  to  concentration.  Fig.  2(b)  shows  the  median  raw  signal  detected  from  each 
nanophosphor  concentration.  Fig.  2(c)  demonstrates  this  method’s  efficacy  in  separating  the 
nanophosphors;  the  correlation  between  the  samples  was  (r  =  -0.98). 
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Fig  2.  (a)  Shown  are  the  Bao.ssYo.iFiiTb’*  and  Bao  ssYo  iFtiEu^*  phosphor  concentrations  with  increasing/decreasing 
concentration  (top-to-bottom)  of  Bao.55Yo.3F2:Tb^^  /  Bao.55Yo.3F2:Eu’'’',  respectively,  (b)  The  raw  signal  detected  for 
each  respective  nanophosphor,  (c)  The  relative  median  recovered  concentration  in  each  ROI  plotted  with  respect  to 

concentration. 


3.2  Animal  phantom  validation 

The  locations  of  the  subcutaneously  injected  RLNPs  are  shown  in  Fig.  3(a).  The  feasibility  of 
this  method  in  a  pre-clinical  subject  with  X-ray  excitation  is  shown  in  Fig.  3(b)-(d).  Tb^’’^- 

doped  RLNPs  embedded  in  the  left  and  right  flanks  are  shown  successfully  recovered  in  Fig. 
3(b).  The  Eu^^-doped  RLNPs  in  the  left  foreleg  and  the  right  flank  are  shown  successfully 
recovered  in  Pig.  3(c).  As  a  control,  inactive  RLNPs  in  the  right  foreleg  did  not  luminesce. 
Pig.  3(d)  shows  the  results  of  the  spectral  unmixing  algorithm,  with  the  Eu^^-doped  RLNPs 
shown  in  red,  and  the  Tb^^-doped  RLNP  shown  in  green.  The  yellow  in  the  multiplexed 
image  indicates  the  presence  of  both  types  of  nanophosphors. 

The  results  of  RLI  with  radiopharmaceutical  excitation  are  shown  in  Pig.  4  (performed 
several  days  prior  but  on  the  same  mouse  as  Pig.  3).  The  PET  image  of  the  four  injected 
regions  of  **F  with  active  and  inactive  particles  is  shown  in  Pig.  4(a);  this  image  is  included 
for  spatial  validation.  All  corners  illuminate  in  the  PET  image  because  of  the  injection  of  ’*P 
with  all  particles.  The  unmixed  Cerenkov  image  is  shown  in  Pig.  4(b).  Note  the  Cerenkov 
signal  in  each  region,  due  to  the  injection  of  ’*P  in  each  region;  these  correspond  to  the 


regions  in  the  PET  image.  In  particular,  note  the  recovery  of  the  Cerenkov  light  in  the  right 
foreleg,  which  contained  the  inactive  control  Bao.55Yo.3F2  nanocrystals  along  with  '*F.  The 
successful  recovery  of  Tb^'^-doped  RFNPs  and  Eu^'^-doped  RFNPs  is  demonstrated  in  Fig.  4c- 
d.  The  spectrally  unmixed  fused  image  is  shown  in  Fig.  4e,  with  Eu^'^-doped  RFNPs  in  red, 
Tb^'^-doped  RFNPs  in  green,  and  Cerenkov  light  in  blue.  Note  that  the  doped  nanoparticles 
did  not  emit  in  the  control  foreleg. 


Fig  3.  X-ray  Luminescence:  (a)  Schematic  of  the  locations  of  each  type  of  RLNP.  The  inactive  particles  are  indicated 
with  the  abbreviation  (Ctrl).  (b,c)  The  unmixed  signal  from  the  Bao  55Yo.3F2:Tb  and  Bao  55Yo.3F2:Eu  particles, 
respectively,  (d)  The  unmixed  multiplexed  image  with  colorbars  for  the  relative  concentrations  of  Bao  55Yo.3F2:Tb  and 

Bao.55Yo..3F2:Eu. 


Fig  4.  Radiopharmaceutical  luminescence:  (a)  Positron  emission  tomography  image  of  the  relevant  RLNPs.  (b)  The 
unmixed  signal  from  Cerenkov  emission.  (c,d)  Unmixed  signal  from  the  Bao.55Yo.3F2:Tb“’*  particles  and 
Bao.55Yo.3F2:Eu^*p^''ticles,  respectively,  (e)  Unmixed  multiplexed  image  with  colorbars  for  the  relative  contributions 
of  Bao.55Yo.3F2:Tb^*,  Bao.55Yo.3F2:Eu’^,  and  Cerenkov  emission. 

4.  Discussion  and  conclusions 

The  methods  outlined  in  this  paper  demonstrate  the  ability  of  this  particle  system  to 
simultaneously  be  used  for  multiple  types  of  RFI;  these  may  then  be  attached  to  distinct 
molecular  targets.  RFNPs  may  become  an  Important  analog  to  the  current  tools  dedicated  to 
pre-clinical  imaging  due  to  their  great  versatility.  In  one  application,  if  injected  with  a 
radiopharmaceutical,  the  nanoparticles  introduced  in  this  study  would  enable  multiple¬ 
targeting;  this  could  prove  to  be  an  Important  method  to  identify  the  kinetics  or  accumulation 
after  drug-injection.  These  RFNPs  may  also  aid  in  the  down-conversion  of  Cerenkov  light 
emitted  from  the  radiopharmaceutical  [14,  15,  16,  17].  In  another  application,  these  RFNPs 
may  be  used  in  conjunction  with  X-ray  luminescence  computed  tomography  [12]  to  recover 
hlgh-resolution  molecular-specific  images.  In  yet  another  application,  this  RFNP  system  may 
aid  in  the  on-line  identification  of  response  to  radiation  therapy.  Through  targeting  of 


molecular  markers  specific  for  tumors,  this  method  could  provide  feedback  during  therapy  as 
to  the  status  of  treatment;  the  penetration  of  high-energy  photons  in  conjunction  with  a  light 
detector  would  ensure  the  ability  to  reach  a  larger  depth  in  tissue. 

In  conclusion,  we  have  demonstrated  a  nanoparticle  platform  for  radioluminescence 
imaging  which  may  be  used  to  image  multiple  nanophosphors  simultaneously  using  both 
radiopharmaceutical  and  X-ray  irradiation.  We  showed  an  inverse  relationship  between 
concentrations  of  the  types  of  RLNPs,  which  showed  successful  separation  of  the  contribution 
from  the  two  emission  signals.  We  then  demonstrated  this  technology  in  a  small-animal 
model. 
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20  CCD  camera  and  a  spectrophotometer;  dose  measurements  were  determined  with  a  calibrated 

21  dosimeter.  Using  these  properties,  in  addition  to  luminescence  efficiency  values  found  in  the  litera- 

22  ture  for  a  similar  phosphor,  minimum  concentration  calculations  are  performed.  Finally,  a  2.5  cm 

23  agar  phantom  with  a  1  cm  diameter  cylindrical  phosphor-filled  inclusion  (diluted  at  10  mg/ml)  is 

24  imaged  to  compare  x-ray  luminescence  contrast  with  x-ray  fluoroscopic  contrast  at  a  superficial 

25  location. 

26  Results:  Dose  to  the  CCD  camera  in  the  chosen  imaging  geometry  was  measured  at  less  than  0.02 

27  cGy/s.  Emitted  light  was  found  to  be  linear  with  dose  (/?^=1)  and  concentration  (R^=l).  Emission 

28  peaks  for  clinical  x-ray  energies  are  less  than  3  nm  full  width  at  half  maximum,  as  expected  from 

29  lanthanide  dopants.  The  minimum  practical  concentration  necessary  to  detect  luminescent  phos- 

30  phors  is  dependent  on  dose;  it  is  estimated  that  subpicomolar  concentrations  are  detectable  at  the 

31  surface  of  the  tissue  with  typical  mammographic  doses,  with  the  minimum  detectable  concentration 

32  increasing  with  depth  and  decreasing  with  dose.  In  a  reflection  geometry,  x-ray  luminescence  had 

33  nearly  a  430-fold  greater  contrast  to  background  than  x-ray  fluoroscopy. 

34  Conclusions:  X-ray  luminescence  has  the  potential  to  be  a  promising  new  modality  for  enabling 

35  molecular  imaging  within  x-ray  scanners.  Although  much  work  needs  to  be  done  to  ensure  bio- 

36  compatibility  of  x-ray  exciting  phosphors,  the  benefits  of  this  modality,  highlighted  in  this  work, 

37  encourage  further  study.  ©  2010  American  Association  of  Physicists  in  Medicine. 

38  [DOT  10.1118/1.3457332] 


39  I.  INTRODUCTION 

AQ:  40  Molecular  imaging  promises  increased  sensitivity  and  speci- 

41  ficity  to  disease  compared  to  traditional  anatomical  imaging 

42  modalities.  The  information  gained  from  molecular  imaging 

43  has  the  potential  to  provide  patient-specific  selection  of 

44  therapy,  improved  prediction  of  outcomes,  and  increased 

45  treatment  efficacy.*  X-ray  radiography  and  computed  tomog- 

46  raphy  (CT)  are  commonly  used  anatomical  imaging  modali- 

47  ties;  however,  although  they  provide  invaluable  information 

48  in  the  clinic,  they  have  been  largely  unsuccessful  for  molecu- 

49  lar  imaging.^  This  deficiency  is  due  to  their  lack  of  sensitiv- 

50  ity  to  low  concentrations  of  contrast  agents;  x-ray  imaging  is 

51  commonly  believed  to  be  many  orders  of  magnitude  less 

52  sensitive  than  opticar  or  radionuclide"*  imaging.  This  poor 


sensitivity  arises  from  the  low  x-ray  stopping  power  of  di- 
luted  contrast  agents,  which  necessitates  high  concentrations  54 
compared  to  other  imaging  modalities.  55 

Phosphors  are  well-established  materials  used  universally  56 
in  cathode  ray  tubes  and  light-emitting  diodes  for  their  abil-  57 
ity  to  emit  light  upon  excitation  by  electrons  or  photons.  58 
Phosphors  consist  of  solid-state  crystals,  which  are  typically  59 
doped  with  transition  metals  or  lanthanide  ions.  These  mate-  60 
rials  form  a  system  optimized  to  capture  higher-energy  radia-  61 
tion  and  emit  downcon verted  energy  as  optical  photons.  In  62 
the  context  of  an  x-ray  scanner,  x-ray  photons  transfer  some  63 
or  all  of  their  energy  to  electrons  in  the  solid-state  crystal  64 
through  Compton  and  photoelectric  interactions.^  These  65 
high-energy  electrons  progressively  lose  energy  through  in-  66 
teractions  with  the  atoms,  leaving  a  track  of  excited  electrons  67 
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behind.  When  the  energy  of  the  excited  electrons  in  the  con- 

69  duction  band  is  reduced  to  approximately  two  to  three  times 

70  the  band  gap,  they  may  migrate  into  the  luminescence  center 

71  of  the  phosphor,  the  dopant,  and  recombine  with  holes  to 

72  emit  light.^  Thus,  this  amplification  process  results  in  effec- 

73  five  quantum  efficiencies  (photons  emitted  divided  by  the 

74  photons  absorbed),  which  can  be  much  greater  than  1.  For 

75  example,  on  the  average,  6000  photons  are  produced  for  each 

76  100  keV  x-ray  photon  absorbed  in  one  particular  gadolinium 

77  oxysulfideiterbium  phosphor.^  The  emitted  light  may  then  be 

78  imaged  by  sensitive  optical  detectors. 

79  This  paper  investigates  the  use  of  nanosized  inorganic 

80  phosphors^  as  potential  biological  contrast  agents  for  medical 

81  imaging  in  a  combined  x-ray/optical  instrument.  The  emis- 

82  sion  from  this  contrast  agent  is  evaluated  to  determine  the 

83  practicality  of  this  new  modality.  The  implications  for  this 

84  x-ray  activated  contrast  agent  are  discussed  with  regard  to  its 

85  potential  to  enable  molecular  imaging  during  fluoroscopy, 

86  x-ray  CT,  or  projection  x-ray  imaging. 

87  II.  METHODS 

88  II.A.  Phosphor  fabrication 

89  Trivalent  europium  (Eu)  or  terbium  (Tb)  activated  gado- 

90  linium  oxysulfides  (GOSs)  or  lanthanum  oxysulfides  (LOSs) 

91  were  synthesized  using  appropriate  rare  earth  nitrates 

92  (99.99%  pure)  with  two  standard  methods:  The  gel-polymer 

93  combustion  process  and  the  combination  capping  process,^ 

94  respectively.  After  preparation,  samples  were  heat  treated  at 

95  500-600  °C  for  1-3  h  to  aid  the  migration  of  the  dopant 

96  into  the  crystal  lattice  stmctures.  Next,  the  powders  were  ball 

97  milled  with  glass  beads  (10  /rm)  in  the  presence  of  the  ap- 

98  propriate  surfactant  for  2-3  h. 

99  II. B.  Spectroscopy  and  imaging  of  phosphor 

100  characteristics 

101  To  facilitate  spectroscopy  and  imaging  for  the  analysis  of 

102  the  properties  of  the  phosphors,  dry  phosphor  was  placed  in 

103  plastic  test  tubes.  For  spectroscopy,  the  distal  end  of  a 

104  10  m-400  yum  optical  fiber  was  placed  in  contact  with  the 

105  side  of  the  test  tube,  while  the  proximal  end  was  attached  to 

106  a  spectrophotometer  (Jaz,  Ocean  Optics,  Dunedin,  FL), 

107  which  was  operated  from  the  x-ray  control  room.  Optical 

108  emission  spectra  across  the  visible  and  near  infrared  (NIR) 

109  range  were  acquired  with  the  SPECTRASUITE  (Ocean  Optics) 

110  software  package.  Imaging  was  performed  with  a 

111  512X512  pixel  backilluminated  CCD  camera  (CCD  tem- 
112perature  maintained  at  -70  °C)  with  a  F1.4  imaging  lens, 

113  exposure  times  varying  from  <1  to  60  s,  and  the  lens  aper- 

114  ture  fully  open.  During  data  acquisition,  the  imaging  camera 

115  was  shielded  with  lead  bricks  and  placed  15-20  cm  outside 

116  the  direct  field  of  radiation  to  protect  from  ionizing  radiation. 

117  The  optics  setup  was  placed  inside  a  light-tight  box  to  elimi- 

118  nate  ambient  room  light.  A  schematic  of  this  imaging  setup  is 

119  shown  in  Fig.  1(a),  and  a  photograph  of  the  setup  during 

120  experimental  imaging  is  shown  in  Fig.  1(b). 


Fig.  1.  (a)  Schematic  of  the  experimental  setup  including  the  CCD  camera, 
the  x-ray  source  and  detector,  and  the  sample,  (b)  Photograph  of  the  imaging 
setup. 

A  cone  beam  computed  tomography  (CBCT)  system 
(Acuity,  Varian  Medical  Systems,  Palo  Alto,  CA)  was  used  to  122 
irradiate  the  sample.  This  system  allows  both  CBCT  and  123 
fluoroscopy  at  various  tube  voltages  between  70  and  125  kV.  124 
This  system  was  operated  in  fluoroscopy  mode  to  enable  125 
continuous  irradiation  from  a  constant  geometry.  126 

II. C.  Concentration  evaluation  127 

Minimum  detectable  concentrations  were  calculated,  as-  128 
suming  a  superficial  location  of  the  phosphor  contrast  agent  129 
(i.e.,  no  signal  loss  due  to  tissue  absorption).  Including  light  130 
detection  losses  L,  the  light  detected  O  is  131 

c^  =  TDcL,  (1)  132 

where  F  is  the  emitted  light  efficiency,  D  is  the  dose,  and  c  is  133 
the  concentration  of  the  phosphor.  From  Ref.  8,  the  emitted  134 
light  efficiency  (in  a  lanthanum  oxysulfideiterbium  phos-  135 
phor)  is  1.39X10'^  optical  photons/ (Gy  mg).  We  assumed  136 
99%  detection  loss  due  to  suboptimal  optical  collection  ge-  137 
ometry.  Signal  below  a  signal-to-noise  ratio  (SNR)  of  10  138 
(assuming  shot-noise  limited  detection)  was  assumed  to  be  139 
too  low  to  detect.  To  calculate  the  molar  concentration,  we  140 
assumed  a  spherical  10  nm  diameter  nanoparticle  consisting  141 
of  hexagonal-stmctured  phosphors  with  lattice  constants  of  142 
fl=4.046  A  and  c=6.951  A  (Ref.  9)  and  density  of  143 
5.5  g/cm^.'°  144 

II. D.  Optical  phantom  fabrication  and  imaging  145 

A  small-animal  sized  tissue-simulating  phantom  was  146 
fabricated  for  this  study.  The  cylindrical  phantom  147 
measured  2.5  cm  in  diameter  X  4.5  cm  in  height,  with  a  148 
1X2.5  cm^  cylindrical  inclusion.  The  phantom  was  made  149 
from  1%  agar  with  homogeneous  optical  properties  using  150 
titanium  oxide  for  scatter  and  India  Ink  for  absorption  using  1 51 
methods  common  to  diffuse  optical  phantoms.**  The  optical  152 
properties  were  determined  by  a  previously  established  153 
system*^  to  be  /u,„=0.0025  mm“*  and  fj./ =0.77  mm“*  at  154 
630.  Micrometer-sized  GOS:Eu  phosphor  particles  were  155 
added  to  the  inclusion  at  a  concentration  of  10  mg/ml,  and  no  156 
phosphor  was  added  to  the  background.  This  phantom  was  157 
imaged  with  an  electron  multiplied  (EM)-CCD  (Hamamatsu  158 
ImageEM  9100-13,  Hamamatsu,  Japan)  with  a  512  159 
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Fig.  2.  (a)  Radiation  dose  on  detector  for  three  different  configurations:  (i)  Direct  irradiation,  (ii)  direct  irradiation  with  lead  shielding  (5  mm  on  top,  50  mm 
on  each  side),  and  (iii)  indirect  irradiation  with  lead  shielding,  (b)  Radiation  events  per  area  per  dose. 


X512  pixel  sensor  cooled  at  -70  °C,  with  a  F1.4  imaging 

161  lens,  with  exposure  of  times  less  of  2  s,  gain  at  half- 

162  maximum,  and  the  lens  aperture  fully  open. 

163  III.  RESULTS 

164  III. A.  CCD  interaction  with  ionizing  radiation 

165  To  ensure  that  the  CCD  camera  would  not  be  damaged 

166  from  the  ionizing  radiation,  an  ionization  chamber  (PTW 

167  Farmer  30010,  PTW  Freiberg  GMBFl,  Germany)  was  placed 

168  in  the  radiation  field  to  determine  the  dose  to  air  in  the  vi- 

169  cinity  of  the  CCD.  Several  locations  along  the  patient  bed 

170  were  measured  to  determine  dose  rate,  as  indicated  in  Fig. 

171  2(a).  Under  direct  radiation,  the  dose  was  nearly 

172  300  yu,Gy/s;  this  dose  rate  was  reduced  by  two  orders  of 

173  magnitude  by  placing  lead  around  the  chamber,  and  further 

174  by  twofold,  to  1.5  /mGy/s,  by  moving  the  chamber  15  cm 

175  out  of  the  radiation  field.  This  rate  deposits  dose  well  below 

176  levels  that  would  likely  damage  the  CCD. 

177  To  investigate  camera  damage  further,  we  investigated  the 

178  lingering  effects  of  radiation  on  the  CCD.  The  interaction  of 

179  an  x-ray  photon  with  the  CCD  camera  appears  in  the  image 

180  as  a  bright  pixel  at  near-maximum  intensity.  These  hot  pixels 

181  appear  similar  to  cosmic  ray  interactions,  which  are  common 

182  with  CCD  cameras.  We  investigated  the  incidence  of  these 

183  events  for  a  typical  setup  using  a  CCD  to  radiation  field 

184  distance  of  15  cm  and  6  mm  of  lead  shielding  above  the 

185  camera  to  protect  from  x-ray  collimator  leakage.  Radiation 

1 86  events  were  determined  by  performing  an  intensity  threshold 

187  on  an  image  acquired  with  the  lens  cap  on.  It  is  clear  from 

188  Fig.  2(b)  that  the  number  of  radiation  events  is  linear  with 

189  dose,  and  thus  there  were  no  lingering  effects  from  the  ra- 

190  diation.  To  improve  image  quality,  denoising  strategies  may 

191  be  employed  utilizing  this  linearity,  such  as  an  automatic 

192  selection  of  a  hot-pixel  threshold,  which  is  dependent  on 

193  camera  dose.  Further  consideration  of  the  damage  limits  of 

194  CCD  cameras  is  given  in  Sec.  IV. 


III.B.  Phosphor  characterization 

III.B.1.  Spectral  emission  196 

A  large  body  of  knowledge  exists  on  phosphors  due  to  197 
over  a  half-century  of  study  of  optimizing  phosphors  for  198 
such  applications  as  light-emitting  diodes,  cathode  ray  tubes,  199 
and  scintillators  for  medical  imaging.  This  work  has  resulted  200 
in  a  library  of  crystals  and  dopants  from  which  one  may  201 
select  an  emission  wavelength  that  is  ideal  for  a  particular  202 
application.*’  We  investigated  the  feasibility  of  GOS  phos-  203 
phors,  which  were  doped  with  either  terbium  (GOS:Tb)  or  204 
europium  (GOS:Eu),  because  of  their  absorption  K-edge  in  205 
the  diagnostic  energy  regime  at  approximately  50  keV.*'^  Fig-  206 
ure  3(a)  shows  the  emission  of  these  phosphors  under  100  207 
kV  x-ray  irradiation.  The  GOS:Tb  phosphor  had  a  maximum  208 
peak  emission  of  545  nm  in  green,  whereas  the  GOS:Eu  209 
phosphor  had  several  peaks  of  longer  wavelengths  in  the  210 
NIR,  including  596,  618,  627,  and  707  nm,  with  an  emission  211 
maximum  at  627  nm.  The  flexibility  enabled  by  modifying  212 
the  dopant  is  of  great  value  for  matching  the  emission  wave-  213 
length  to  a  particular  application,  such  as  the  absorption  peak  214 
of  a  phototherapeutic  brng,*"*  or  the  tissue  absorption  mini-  215 
mum  for  optical  imaging  in  deep  tissue.'^  216 


III.B.2.  Light  output  vs  dose  217 

To  determine  the  linearity  of  light  output  from  phosphor,  218 
GOS:Tb  was  dispersed  in  a  cuvette  containing  1%  agar  at  a  219 
concentration  of  10  mg/ml.  The  phosphor  solution  was  220 
placed  at  the  same  source-target  distance  as  the  ionization  221 
chamber.  The  x-ray  system  was  operated  in  fluoroscopy  222 
mode  and  the  tube  voltage  was  set  to  lOOkV.  The  dose  was  223 
linearly  increased  by  two  methods.  First,  the  tube  current  224 
was  increased  from  5  to  20  mA  with  a  constant  tube  voltage.  225 
Second,  the  tube  voltage  and  current  remained  constant,  and  226 
dose  was  linearly  increased  by  adjusting  the  irradiation  time.  227 
Images  were  acquired  with  a  CCD  camera,  with  the  exposure  228 
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Fig.  3.  (a)  Emission  spectra  of  GOS  phosphors  under  x-ray  excitation.  GOS:Tb,  with  an  emission  peak  at  545  nm,  is  labeled  with  asterisks  (*),  and  GOS:Eu, 
with  emission  peaks  at  596,  618,  627,  and  707  nm.  (b)  X-ray  induced  GOSiTb  Photon  output  vs  dose  with  dose  varied  by  increasing  tube  current,  (c)  Photon 
output  vs  dose  with  dose  varied  by  increasing  irradiation  time.  For  both  methods,  photon  output  is  linear  with  dose  shown  in  terms  of  current 
(r=  1 ,  p<  0.001 ,  R^=  1).  (d)  Photon  output  vs  concentration  for  the  GOS:Eu  phosphor.  Concentrations  were  measured  via  microcuvette  and  dispersed  in  1% 
agar.  Photon  output  is  linear  with  concentration  (r=0.99,  p=1.9X  10“^). 


time  optimized.  All  measurements  were  denoised  for  hot 

230  spots  and  identical  regions  of  interest  (ROIs)  were  selected 

231  for  analysis. 

AQ:  232  Both  current/dose  and  time/dose  linearity  were  confirmed 

233  by  mnning  a  linear  correlation  analysis  (correlation  coeffi- 

234  cient  of  1.0,  p<0.001).  Figure  3(b)  shows  the  linearity  in 

235  dose  using  the  first  method,  which  showed  significant  linear- 

236  ity  (correlation  coefficient  of  1.0,  p<  0.001),  while  Fig.  3(c) 

237  confirms  linearity  with  the  second  method  (correlation  coef- 

238  ficient  of  1.0,  p<0.001). 

239  III.B.3.  Light  output  vs  concentration 

240  To  assess  the  light  output  due  to  various  phosphor  con- 

241  centrations,  dilutions  of  GOS:Eu  phosphors  from  5  /rg/ml 

242  to  10  mg/ml  were  dispersed  in  1%  agar.  Cuvettes  were 

243  placed  in  a  50  kV,  30  mA  x-ray  source  and  imaged  with  a 

244  CCD  camera.  ROIs  were  selected  to  contain  similar  areas  of 

245  the  cuvettes,  and  the  signal  was  normalized  according  to  ex- 

246  posure  time.  Hot  spots  were  removed  from  the  images  prior 

247  to  analysis. 

248  The  results  shown  in  Fig.  3(d)  demonstrate  a  strong  lin- 

249  earity  (linear  correlation  coefficient  of  0.99,  p<  0.001)  with 

250  a  slope  of  1.09.  The  slight  departure  of  the  slope  from  unity 

251  is  most  likely  due  to  errors  in  selecting  identical  ROIs  for  the 


different  cuvettes  and  may  have  resulted  in  the  inclusion  of 
the  cuvette  wall  in  the  ROI,  which  exhibited  some  light  pip-  253 
ing.  254 

III.C.  Minimum  detectable  concentrations  255 

We  calculated  the  minimum  detectable  concentration  ac-  256 
cording  to  the  methodology  outlined  in  Sec.  II C,  for  doses  257 
ranging  from  1  cCy  (less  than  the  typical  mammographic  258 
dose)  to  20  Cy  (a  typical  dose  delivered  in  single-dose  intra-  259 
operative  radiation  therapy).  In  addition,  we  calculated  the  260 
minimum  detectable  concentration  for  several  phosphor  effi-  261 
ciencies,  scaled  according  to  that  reported  by  Kandarakis  et  262 
al.^  (i.e.,  100%  is  equivalent  to  the  efficiency  reported).  Fig-  263 
ure  4  shows  the  minimum  detectable  concentrations  for  these  264 
scenarios.  According  to  these  calculations,  picomolar  (ng/ml)  265 
concentrations  are  detectable  (SNR  of  10)  with  a  LOS  phos-  266 
phor  for  mammographiclike  dose,  while  therapy  doses  allow  267 
femtomolar  (pg/ml)  concentrations  to  be  detected.  268 

III.D.  Contrast  comparison  between  x-ray/optical  and  269 
fluoroscopy  270 

The  recovered  contrast  between  an  inclusion  with  phos-  271 
phor  and  a  background  without  phosphor  was  investigated  to  272 
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Fig.  4.  Minimum  estimated  phosphor  concentrations  calculated  from  the  literature  and  our  own  data.  Lines  represent  phosphor  luminescence  efficiencies;  the 
100%  line  is  based  on  efficiency  data  from  a  La202S  :Tb  phosphor  reported  by  Kandarakis  et  al.  (a)  Minimum  weight  concentration  (in  mg/ml)  vs  dose  (Gy), 
(b)  Minimum  molar  concentration  (in  M)  vs  dose  (Gy). 


compare  the  contrast  differences  between  x-ray  fluoroscopy 

274  and  x-ray/optical  luminescence  imaging  in  the  small-animal 

275  imaging  phantom  described  in  Sec.  II D.  The  phantom  was 

276  imaged  during  fluoroscopy  operation  with  the  tube  voltage  at 

277  100  kV  and  the  tube  current  at  10  mA.  The  x-ray  source  was 

278  placed  above  the  phantom  to  evenly  irradiate  the  volume. 

279  The  CCD  camera  was  placed  within  15  cm  of  the  phantom 

280  and  oriented  orthogonal  to  the  direction  of  irradiation.  The 

281  fluoroscopy  image  was  taken  simultaneously.  It  is  important 

282  to  note  that  gadolinium  is  a  strong  absorber  of  100  keV  x-ray 

283  energy,  with  a  mass  attenuation  coefficient  of  3.109  cm^/g 

284  (compared  to  common  x-ray  contrast  agents  such  as  barium 

285  at  2.196  cm^/g  and  iodine  at  1.942  cm^/g).  Since  gado- 

286  linium  has  a  higher  mass  attenuation  coefficient  for  x-ray 

287  photons  than  water  (mass  attenuation  coefficient  of 

288  0.1707  cm^/g),  it  should  exhibit  slight  contrast. 

289  The  images  from  the  phantom  are  shown  in  Fig.  5.  Figure 

290  5(a)  shows  a  white  light  image  taken  by  the  CCD  camera 

291  with  background  illumination  from  the  room  lights.  The  cor- 

292  responding  fluoroscopy  image  is  shown  in  Fig.  5(b).  The 

293  phosphor  inclusion  is  indicated  by  the  red  arrow  and  shows 

294  slight  increased  x-ray  absorption  (the  smaller  dark  circle  is  a 

295  bolt  hole  in  the  aluminum  optical  table).  Figure  5(c)  shows  a 

296  raw  optical  image  taken  of  the  x-ray  luminescent  phosphor. 

297  This  image  is  overlaid  on  the  white  light  image  in  Fig.  5(d). 

298  The  contrast  between  the  inclusion  and  the  background  is 

299  very  slight,  at  0.6%  for  the  fluoroscopy  image,  while  it  is 

300  over  260%  for  the  luminescence  image.  In  addition,  the 

301  signal-to-noise  ratio  for  the  phosphor  emission  was  23  vs  2.4 

302  for  the  fluoroscopy  image. 

303  IV.  DISCUSSION 

304  We  found  that  dose  distributed  to  the  shielded  camera  was 

305  measured  at  less  than  3  yuGy/s  when  the  camera  was  posi- 

306  tioned  at  the  isocenter  and  the  x-ray  tube  voltage  was  100 

307  keV  with  the  tube  current  at  20  mA.  Although  this  dose  is 

308  low,  estimating  a  damage  threshold  is  difficult  for  CCDs  be- 

309  cause  damage  is  design/manufacturer-speciflc  and  is  depen- 
31  0  dent  on  environmental  conditions  (for  a  more  thorough  over- 
311  view,  see  Ref.  16).  It  is  well  recognized  that  the  largest 


radiation  threat  to  the  operation  of  a  CCD  is  the  bombard- 
ment  by  highly  energetic  heavy  particles,  such  as  protons  and  313 
neutrons.  These  particles  contribute  most  to  CCD  damage  314 
through  impact  displacements  of  silicon  atoms  which  create  315 
semipermanent  energy  traps.  These  traps  create  energy  levels  316 
which  can  increase  Johnson  noise  (via  promoting  valence  317 
band  electrons  to  the  conduction  band),  create  spurious  noise  318 
when  trapped  electrons  are  released,  and  alter  the  operation  319 
of  transistor  gates  by  altering  their  flat-band  voltage  (for  a  320 
more  thorough  review,  see  Refs.  17  and  18).  Although  the  321 
probability  of  creating  protons  and  neutrons  is  extremely  low  322 
at  diagnostic  x-ray  energies  studied  here,  it  is  relevant  for  323 
therapeutic  energies  in  the  MV  range.  Cumulative  doses  are  324 
also  important  because  of  the  increased  probability  for  a  325 
high-energy  photon  interaction.  It  has  been  reported  that  total  326 


Fig.  5.  Agar  imaging  phantom  with  embedded  phosphors  and  tissue  emu¬ 
lating  optical  properties,  (a)  White-light  optical  image,  (b)  Projection  fluo-  AQ: 

roscopy  image  (note  the  distinction  between  the  phosphor  inclusion —  ^12 

indicated  by  the  red  arrow — around  300  units,  compared  to  the  black  circle 
caused  by  a  screw  hole  in  the  optical  table  supporting  the  phantom),  (c) 

Optical  emission  from  the  phantom,  (d)  Overlay  of  the  white  light  image  (a) 
and  the  light  emission  (c). 
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doses  above  about  25  and  100  Gy  are  considered  thresholds 

328  for  increased  noise  and  permanent  damage,  respectively*® 

329  (though  again,  these  numbers  are  highly  camera-dependent). 

330  Thus,  it  is  important  to  keep  radiation  dose  on  the  CCD  as 

331  low  as  possible  to  minimize  cumulative  dose  and  increase 

332  the  life  of  the  camera.  The  dose  deposited  in  this  experiment 

333  of  1-3  /xGy/s  is  minimal,  and  damage  due  to  radiation 

334  should  be  insignificant  in  the  lifetime  of  the  camera. 

335  The  lanthanide  dopants  examined  in  this  study  are  known 

336  to  have  extremely  sharp  peaks;  for  example,  we  measured 

337  the  full  width  at  half  maximum  (FWHM)  of  the  GOS:Tb 

338  sample  at  3  nm,  while  the  GOS:Eu  sample  had  a  FWHM  of 

339  2  nm  for  the  largest  peak  [both  shown  in  Fig.  3(a)].  This  is  in 

340  contrast  to  common  fluorophores  like  indocyanine  green, 

341  which  have  often  have  emission  peak  FWHM  of  about  60 

342  nm.  This  enables  the  possibility  of  multiplexing  contrast 

343  agents  with  very  little  cross-talk,  allowing  the  simultaneous 

344  measurement  of  several  markers  for  disease. 

345  A  major  concern  of  these  x-ray  excitable  phosphor  con- 

346  trast  agents  is  that  they  require  ionizing  dose  to  activate. 

347  Thus,  lower  concentrations  of  phosphors  will  necessitate 

348  higher  doses.  We  analyzed  the  practical  detection  limit  of 

349  phosphors  using  the  knowledge  of  concentration  and  dose 

350  linearity  and  luminescence  efficiencies  found  in  the  literature 

351  for  similar  nanophosphors.  For  detection  at  the  surface  of 

352  tissue,  x-ray  mammographic  dose  should  be  sufficient  to  al- 

353  low  the  detection  of  picomolar  (ng/ml)  concentrations  of 

354  phosphor.  This  finding  is  corroborated  with  our  experimental 

355  results  shown  in  Fig.  3.  If  the  data  in  Fig.  3(d)  are  extrapo- 

356  lated  to  10  counts/pixel,  adequate  signal-to-noise  ratio  is 

357  achieved  using  an  FM-CCD  camera  (with  dark  noise  of  0.05 

358  counts/pixel  for  this  acquisition).  Placing  the  camera  within  a 

359  few  cm  (compared  to  20  cm  in  this  experiment)  would  result 

360  in  the  detection  of  approximately  two  orders  of  magnitude 

361  more  light.  Further,  if  the  dose  is  increased  by  an  order  of 

362  magnitude  [from  the  4  mGy  dose  administered  for  Fig.  3(d)] 

363  to  mammography  levels,  and  FM  gain  is  increased  to  the 

364  maximum,  concentrations  of  ng/ml  (picomolar)  may  be  real- 

365  ized. 

366  At  deeper  depths,  however,  light  attenuates  rapidly;  for 

367  example,  in  breast  tissue,  a  depth  of  3  cm  would  attenuate 

20 

368  detected  light  by  approximately  three  orders  of  magnitude. 

369  This  would  effectively  decrease  phosphor  sensitivity  during 

370  mammographic  dose  to  nanomolar  concentrations.  In  a  ra- 

371  diation  oncology  setting,  however,  high  doses  are  delivered 

372  to  treat  disease  (such  as  the  20  Gy  single-fraction  therapy, 

373  which  is  used  during  intraoperative  radiation  therapy^*).  This 

374  technique  could  aid  a  surgeon  and  radiation  oncologist  to 

375  identify  the  distribution  of  disease  around  a  tumor  margin, 

376  such  as  during  breast  lumpectomy.  In  this  case,  the  patient 

377  would  be  injected  with  a  phosphor  contrast  agent  before  or 

378  during  surgery,  depending  on  the  pharmacokinetics  of  the 

379  tracer.  The  tracer  could  be  imaged  during  the  first  1%  or  10% 

380  (e.g.,  0.2  or  2  Gy)  of  the  radiotherapy  treatment  dose,  which 

381  would  provide  the  clinicians  with  more  confidence  about  the 

382  treatment  dose  or  volume,  or  enable  an  adjustment  to  the 

383  dose  distribution. 


The  other  area  of  concern  with  these  particles  is  biologi- 
cal  compatibility.  This  is  an  issue  for  all  nanoparticle  sys-  385 
terns,  and  much  work  is  being  done  to  develop  strategies  to  386 
ensure  stability  and  compatibility.^^  In  fact,  multiple  groups  387 
have  successfully  used  upconversion  phosphors  in  biological  388 
small-animal  models. The  increased  interest  in  phosphors  389 
should  aid  in  the  rapid  advancement  in  biocompatibility,  390 
which  will  aid  this  modality.  391 

We  found  that  the  contrast  to  background  ratio  was  over  392 
2.5  orders  of  magnitude  higher  for  optical  detection  of  the  393 
luminescent  phosphors  compared  to  x-ray  fluoroscopy.  The  394 
actual  contrast  amplification  is  much  higher  since  the  optical  395 
photons  emitted  from  the  inclusion  exhibit  extensive  scatter  396 
in  the  background  and  subsequently  are  detected  by  the  397 
CCD.  These  scatter  effects  would  be  greatly  minimized  via  398 
modeling  of  the  light  propagation.  In  comparison,  x-ray  pho-  399 
ton  scatter  is  relatively  nonexistent  so  that  the  signal  contri-  400 
bution  from  the  background  originating  from  the  inclusion  is  401 
negligible.  Considering  the  photon  scatter  should  greatly  im-  402 
prove  the  contrast  of  these  phosphors  for  optical  detection  403 
compared  to  fluoroscopy.  404 

Although  this  work  demonstrated  the  potential  of  x-ray  405 
luminescence  imaging  for  imaging  a  superficial  object,  im-  406 
aging  of  lesions  centimeter  deep  should  be  possible,  with  407 
contrast-resolution  limitations  depending  on  tissues’  proper-  408 
ties,  concentration,  and  nonspecific  uptake.  The  development  409 
of  deep-tissue  x-ray  luminescence  imaging  will  require  the  410 
incorporation  of  optical  tomographic  models.  With  x-ray  lu-  411 
minescent  imaging,  the  x-ray  source  must  be  modeled  in  412 
tissue  to  give  an  accurate  description  of  dose.  There  are  413 
many  sophisticated  tools  to  model  dose,  such  as  Monte  Carlo  414 
or  analytical  models,  which  have  been  shown  to  be  accurate  415 
[e.g.,  within  4%  (Ref.  25)]  in  biological  tissues.  Concurrent  416 
x-ray  stmctural  imaging  will  further  improve  these  calcula-  417 
tions.  After  dose  distribution  is  calculated,  tomographic  im-  418 
aging  may  be  performed  with  a  reconstruction  model  that  419 
uses  a  model  of  the  light  propagation  in  tissue  to  minimize  420 
the  difference  between  calculated  and  optical  measurements.  421 
This  is  very  similar  to  the  fluorescence  molecular  imaging  422 
problem.^®  Once  again,  the  knowledge  of  anatomical  infer-  423 
mation  will  aid  the  optical  reconstmction  problem  by  provid-  424 
ing  structural  detail  which  may  be  used  to  improve  optical  425 
modeling and  reconstruction.^*  426 

The  joint  use  of  x-ray  activated  phosphors  for  molecular  427 
imaging  offers  several  advantages  to  x-ray  imaging  and  to  428 
all-optical  fluorescent  imaging.  For  x-ray  imaging,  contrast  429 
agent  imaging  is  currently  limited  to  high  concentrations  of  430 
nonspecific  iodine  or  barium  sulfate.  Optimal  concentrations  431 
for  these  contrasts  have  been  reported  around  300-500  432 
mg/ml. These  high  concentrations  are  impractical  for  im-  433 
aging  biological  targets.^  The  ability  to  image  cellular  targets  434 
would  be  a  great  benefit  to  x-ray  imaging,  which,  despite  435 
being  the  most  prominent  modality  in  use  in  the  clinic  today,  436 
is  generally  limited  to  imaging  stmctural  anatomy.  The  use  437 
of  phosphors  combined  with  the  sensitivity  of  optical  imag-  438 
ing  allows  lower,  more  biologically  feasible  concentrations  439 
of  contrast  agents  than  is  currently  available  with  x-ray  im-  440 
aging  alone.  441 
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The  use  of  x-ray  activated  phosphors  offers  three  unique 

443  advantages  to  all-optical  approaches.  First,  this  dual- 

444  modality  instrument  offers  inherent  spatial  coregistration  be- 

445  tween  anatomical  features  and  optical  contrasts.  This  regis- 

446  tration  is  critical  for  imaging  functional  pathology  in  deep 

447  tissue  [hence,  the  need  for  positron  emission  tomography 

448  (PET)/CT  imaging  systems].^**  Next,  the  use  of  x-ray  excita- 

449  tion  eliminates  the  optical  autofluorescence  issue  in  optical 

450  imaging.  Since  the  x-ray  excitation  spectrum  is  undetectable 

451  with  photo-optical  detectors,  autofluorescence  is  avoided, 

452  which  potentially  reduces  the  detection  limit  for  low  concen- 

453  trations.  Finally,  this  technique  is  also  expected  to  have  in- 

454  creased  depth  performance  over  optical  imaging,  because  of 

455  the  high  penetration  of  x-ray  photons  in  tissues.  X-ray  pho- 

456  tons  have  nearly  two  orders  of  magnitude  lower  effective 

457  attenuation  coefficient  compared  to  optical  photons;  this  op- 

458  portunity  offers  the  potential  to  use  clinically  available  in- 

459  strumentation  as  an  external  source. 

460  V.  CONCLUSIONS 

461  This  study  focused  on  the  instmmentation  and  material 

462  feasibility  of  inorganic  downconversion  phosphors  toward 

463  the  realization  of  x-ray  molecular  imaging.  Significant  recent 

464  advances  in  PET,"^  optical  imaging,"  magnetic  resonance 

465  imaging,*^  and  to  a  lesser  extent,  single  positron  emission 

466  computed  tomography,  and  ultrasound  have  invigorated  the 

467  search  for  disease-specific  protein  receptors  that  may  be  tar- 

468  geted  with  imaging  agents.  This  approach  has  already  been 

469  applied  to  numerous  pathologies  to  identify  atherosclerosis 

470  and  thrombosis,  to  determine  treatment  efficacy  via  apop- 

471  tosis  markers, to  identify  cancer,  and  to  monitor  cellular 

472  activity.  The  incorporation  of  these  markers  into  x-ray  imag- 

473  ing  may  have  significant  impact  on  medical  imaging. 

474  In  this  work,  we  demonstrate,  for  the  first  time  to  our 

475  knowledge,  the  feasibility  of  using  inorganic  phosphors  to 

476  enable  optical  detection  under  x-ray  irradiation,  which  may 

477  enable  x-ray  molecular  imaging.  We  first  investigated  the 

478  practical  feasibility  of  operating  a  CCD  within  an  x-ray  ex- 

479  citation  field  at  clinically  relevant  energies,  taking  into  con- 

480  sideration  noise  and  potential  damage.  We  found  that  the 

481  dose  distribution  to  air  was  sufficiently  low  to  prevent  dam- 

482  age  during  operation.  Additionally,  the  noise  on  the  CCD  due 

483  to  incoming  high-energy  irradiation  is  manageable.  We  dem- 

484  onstrated  the  ability  of  several  phosphors  to  emit  light  in  the 

485  optical  regime  under  x-ray  excitation.  These  phosphors 

486  should  be  effective  for  tuning  light  output  to  a  specific  ap- 

487  plication.  We  found  that  the  light  output  was  linearly  propor- 

488  tional  to  both  dose  and  concentration.  Future  work  will  focus 

489  on  quantitative  imaging.  We  calculated  minimum  detectable 

490  concentrations  based  on  these  data  and  values  found  in  lit- 

491  erature;  these  concentrations  are  sufficient  for  certain  bio- 

492  logical  imaging  applications.  Finally,  we  demonstrated  the 

493  potential  of  inorganic  phosphors  to  image  lower  concentra- 

494  tions  than  is  possible  with  x  ray  alone.  We  found  a  430  times 

495  improvement  in  contrast  recovery  for  optical  detection  com- 

496  pared  to  fluoroscopic  detection.  This  improvement  is  ex- 

497  pected  to  be  greater  with  modeling  of  photon  propagation 


and  imaging.  We  envision  hybrid  x-ray/optical  imaging  may 
have  significant  application  in  the  detection  and  diagnosis  of  499 
disease,  especially  during  image-guided  intervention.  500 
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Cerenkov  luminescence  imaging  (CLI)  is  an  emerging  new 
molecular  imaging  modality  that  is  relatively  inexpensive,  4s 
easy  to  use,  and  has  high  throughput.  CLI  can  image 
clinically  available  PET  and  SPECT  probes  using  optical 
instrumentation.  Cerenkov  luminescence  endoscopy  (CLE)  is 
one  of  the  most  intriguing  applications  that  promise  potential 
clinical  translation.  We  developed  a  prototype  customized 
fiberscopic  Cerenkov  imaging  system  to  investigate  the 
potential  in  guiding  minimally  invasive  surgical  resection. 
Methods:  All  experiments  were  performed  in  a  dark  cham¬ 
ber.  Cerenkov  luminescence  from  '®F-FDG  samples  contain¬ 
ing  decaying  radioactivity  was  transmitted  through  an  optical 
fiber  bundle  and  imaged  by  an  intensified  charge-coupled 
device  camera.  Phantoms  filled  with  ''®F-FDG  were  used  to 
assess  the  imaging  spatial  resolution.  Finally,  mice  bearing 
subcutaneous  C6  glioma  cells  were  injected  intravenously 
with  '^F-FDG  to  determine  the  feasibility  of  in  vivo  imaging. 
The  tumor  tissues  were  exposed,  and  CLI  was  performed  on 
the  mouse  before  and  after  surgical  removal  of  the  tumor 
using  the  fiber-based  imaging  system  and  compared  with 
a  commercial  optical  imaging  system.  Results:  The  sensitiv¬ 
ity  of  this  particular  setup  was  approximately  45  kBq  (1.21 
(jLCi)/300  [jlL.  The  3  smallest  sets  of  cylindric  holes  in  a  com¬ 
mercial  SPECT  phantom  were  identifiable  via  this  system, 
demonstrating  that  the  system  has  a  resolution  better  than 
1.2  mm.  Finally,  the  in  vivo  tumor  imaging  study  demon¬ 
strated  the  feasibility  of  using  CLI  to  guide  the  resection  of 
tumor  tissues.  Conclusion:  This  proof-of-concept  study  ex¬ 
plored  the  feasibility  of  using  fiber-based  CLE  for  the  detection 
of  tumor  tissue  in  vivo  for  guided  surgery.  With  further  improve¬ 
ments  of  the  imaging  sensitivity  and  spatial  resolution  of  the  cur¬ 
rent  system,  CLE  may  have  a  significant  application  in  the  clinical 
setting  in  the  near  future. 
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V_^erenkov  luminescence  imaging  (CLI)  has  recently  [aqi] 
attracted  increasing  interest  in  the  field  of  molecular  im¬ 
aging  (1,2).  CLI  is  a  new  optical  imaging  modality  in 
which  images  are  obtained  by  monitoring  the  Cerenkov 
photons  emitted  from  highly  energetic  moving  charged 
particles  ((3^  or  p“).  Conventional  nuclear  imaging 
methods,  such  as  PET  and  SPECT,  are  the  most  widely 
used  clinical  molecular  imaging  techniques.  However, 
these  modalities  usually  suffer  from  high  cost,  limited 
availability,  relatively  low  spatial  resolution,  and  low 
throughput  (3).  As  a  technique  that  bridges  optical  imaging 
and  radionuclide  imaging,  CLI  has  shown  many  advantages 
such  as  high  sensitivity,  high  resolution,  low  cost,  wide  avail¬ 
ability,  relatively  high  throughput,  and  commercially  avail¬ 
able  radionuclide  probes  already  approved  by  the  Food  and 
Drug  Administration. 

Since  its  discovery  in  2009,  CLI  has  quickly  become 
a  practical  molecular  imaging  technique,  and  many  new 
applications  of  CLI  in  preclinical  research  continue  to 
emerge  (1,2).  Several  research  groups  have  demonstrated 
that  CLI  can  be  a  powerful  tool  for  tumor  imaging  using 
radionuclide  probes  such  as  '^F-FDG.  Important  validation 
studies  have  also  been  performed,  and  it  was  reported  that 
there  is  a  good  linear  correlation  between  the  tumor  uptake 
quantified  by  PET  and  tumor  CLI  signals  in  subcutaneous 
xenograft  models  (4-7).  CLI  has  also  been  found  to 
be  useful  for  imaging  a  and  pure  |3“  emitters  such  as  ""Y 
and  ^^^Ac,  which  are  used  for  cancer  treatment  (4).  More¬ 
over,  CLI  can  be  applied  in  the  monitoring  of  reporter  gene 
expression.  The  herpes  simplex  virus  type  1  thymidine  ki¬ 
nase  and  the  sodium  iodide  symporter  reporter  genes  were 
recently  successfully  imaged  by  CLI  in  conjunction  with 
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appropriate  radioactive  reporter  probes  (4,8,9).  Cerenkov 
photons  have  also  been  used  to  illuminate  fluorescent 
dyes  and  nanoparticles  for  in  vivo  imaging  (10-13). 
A  recent  study  further  showed  that  CLI  could  be  used 
for  monitoring  tumor  surgery  in  xenograft  mouse  models, 
highlighting  the  translational  potential  of  the  modality 
(14).  CLI  has  also  been  used  to  measure  the  radiochem¬ 
ical  purity  of  a  radiolabeled  compound  and  to  image 
plant  physiology  (15).  In  addition  to  these  applications, 
the  recent  advancement  in  CLI  tomography  makes  this 
novel  imaging  modality  even  more  powerful  and  prom¬ 
ising  for  biomedical  research  including  diagnostic  imag¬ 
ing  and  therapeutic  monitoring  (16-19). 

Because  it  can  image  clinically  available  radiotracers, 
CLI  has  the  potential  to  be  rapidly  translated  into  clinical 
applications  (4,18).  However,  all  the  studies  described 
above  have  involved  the  use  of  conventional  small-animal 
optical  imaging  systems,  which  are  not  compatible  with 
routine  clinical  practice.  In  this  study,  we  have  built  the 
first,  to  our  knowledge,  prototype  system  that  is  amenable 
to  Cerenkov  luminescence  endoscopy  (CLE)  in  the 
clinic.  This  system  comprises  an  optical  fiber  bundle 
and  an  intensified  charge-coupled  device  (CCD)  camera. 
The  small-diameter  flexible  endoscope  is  designed  for  min¬ 
imally  invasive  monitoring  of  living  tissues  and  organs, 
using  a  sensitive  camera  to  produce  images  of  the  radionu¬ 
clide  probe.  It  is  also  expected  that  hollow  organs  (e.g., 
bladder  and  lung)  or  insufflated  cavities  (e.g.,  esophagus 
and  colon)  inside  the  body  can  provide  natural  and  ana¬ 
tomic  dark  chambers  for  CLE.  By  overlaying  conventional 
bright-field  images  commonly  obtained  from  an  endoscope 
with  the  Cerenkov  luminescence  images  generated  by  ra¬ 
dionuclide  probes,  CLE  could  be  used  to  identify  diseased 
tissues  for  diagnostic  purposes  and  real-time  monitoring  of 
endoscopic  surgery. 

Therefore,  we  investigated  the  feasibility  of  CLE  for 
guiding  cancer  surgical  resection  by  performing  sensi¬ 
tivity,  spatial-resolution,  and  proof-of-concept  experi¬ 
ments.  The  characteristic  sensitivity  of  this  system  was 
evaluated  using  the  natural  decay  of  '^E-EDG.  Imaging 
resolution  was  determined  using  a  phantom  with  cylindric 
holes.  Einally,  imaging  studies  were  performed  to  dem¬ 
onstrate  surgical  guidance  in  a  small-animal  tumor  mo¬ 
del.  A  well-known  commercial  optical  imaging  system 
was  used  to  provide  a  comparison.  The  data  presented 
herein  outline  the  instrumentation  and  methodology;  fur¬ 
ther  steps-needed  to  translate  this  modality  into  the  clinic 
are  also  discussed. 

MATERIALS  AND  METHODS 

*®F-FDG  was  produced  by  the  Radiochemistry  Facility  at  Stan¬ 
ford  University.  The  rat  glioma  cell  line  C6  was  obtained  from 
American  Type  Culture  Collection.  Female  athymic  nude  mice 
(nu/nu)  (n  =  5)  obtained  from  Charles  River  Laboratories,  Inc., 
were  4-6  wk  old.  A  CRC-15R  PET  dose  calibrator  (Capintec  Inc.) 
was  used  for  all  radioactivity  measurements. 


Tumor  Model 

All  animal  studies  were  performed  in  compliance  with  federal 
and  local  institutional  guidelines  for  the  conduct  of  animal 
experimentation.  C6  cells  were  cultured  in  Dulbecco  modified 
Eagle  medium  supplemented  with  10%  fetal  bovine  serum  and  1% 
penicillin/streptomycin  (Invitrogen  Life  Technologies).  The  cell 
line  was  maintained  in  a  humidified  atmosphere  of  5%  CO2  at 
37°C,  with  the  medium  changed  every  other  day.  A  75%  confluent 
monolayer  was  detached  with  trypsin  and  dissociated  into  a  sin¬ 
gle-cell  suspension  for  further  cell  culture.  Approximately  1  x  10® 
C6  cells  suspended  in  phosphate-buffered  saline  (0.1  M,  pH  7.2; 
Invitrogen)  were  implanted  subcutaneously  in  the  left  legs  of  nude 
mice.  Tumors  were  allowed  to  grow  to  150-200  mm^  (2-3  wk), 
and  the  tumor-bearing  mice  were  imaged  in  vivo. 

CLI 

Validation  CLI  was  performed  with  an  IVIS  Spectrum  system 
(Caliper  Life  Science).  For  all  in  vivo  studies,  radionuclides  were 
diluted  in  phosphate-buffered  saline.  Animals  were  placed  in 
a  light-tight  chamber  under  isoflurane  anesthesia.  Each  acquisition 
took  3  min  for  all  studies  without  filters.  Images  were  acquired  and 
analyzed  using  Living  Image  3.0  software  (Caliper  life  sciences). 
The  mice  were  kept  fasting  overnight  before  *®F-FDC  imaging  and 
anesthetized  by  inhalation  of  2%  isoflurane  during  the  experiment. 


FIGURE  1.  Suggested  application  fef  fiber-based  system  for  en¬ 
doscopic  and  laparoscopic  CLI. 


Fiber-Based  CLE  System 

The  fiber-based  CLE  system  used  an  imaging  optical  fiber 
coupled  to  a  highly  sensitive  intensified  CCD  camera.  The 
application  of  this  system  for  surgical  oncology  guidance  is 
detailed  in  Figure  1.  Specifically,  a  microimaging  lens  (Cinegon,  [Fig.  1] 
F/1.4,  12-mm  focal  length;  Schneider)  was  coupled  at  the  distal 
end  of  an  optical  imaging  fiber  bundle  that  was  108  mm  long,  with 
a  5  X  6.7  mm  active  area.  The  pixels  of  the  fiber  bundle  are  made 
of  discrete  10-p.m  fibers  (IC154;  Schott).  At  the  proximal  end, 
a  Pentax  (F/1.4)  lens  provided  relay  optics  to  the  camera  (Supple¬ 
mental  Fig.  lA;  supplemental  materials  are  available  online  only 
at  http://jnm.snmjournals.org).  The  camera,  an  image-intensified 
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CCD  (Turbo  640-Z;  Stanford  Photonics  Inc.),  had  640  x  480 
pixels  and  single-photon  imaging  capability  (Supplemental  Fig. 
IB).  To  minimize  background  light  and  to  emulate  a  light-tight 
anatomic  cavity,  all  images  were  taken  in  a  dark  box.  In  addition, 
images  were  postprocessed  with  a  software  thresholding  tech¬ 
nique  that  removed  pixel  values  below  a  manufacturer-recom¬ 
mended  value;  this  optimal  value  was  chosen  to  remove  low- 
intensity  pixels  that  corresponded  to  thermal  and  read  noise  on 
the  CCD.  These  steps  reduced  the  noise  to  approximately  20 
counts  per  second. 

Imaging  Using  Fiber-Based  System 

Images  were  acquired  at  a  high  frame-rate  of  120  Hz  to  allow 
for  maximal  reduction  in  noise.  The  noise-reduction  step  de¬ 
scribed  above  was  performed  for  each  frame.  These  frames  were 
then  accumulated  for  5  min  to  form  the  raw  data.  The  raw  data 
were  image-processed  offline  using  an  algorithm  to  remove 
cosmic  and  stray  y-events.  This  algorithm  removed  noisy  pixels 
by  thresholding  the  gradient  around  each  pixel  in  its  local  area.  An 
additional  denoising  step  removed  spurious  bright  pixels  that 
appeared  in  sequential  exposures. 

The  fiber-based  CLE  system  was  characterized  for  sensitivity 
by  imaging  2  adjacent  wells  of  a  black  96-well  plate  (300  p.L  per 
well).  A  glycerol  and  water  mixture  combined  with  3.7  MBq  (100 
p-Ci)  of  ^*F-FDG  was  in  1  well,  and  a  glycerol  and  water  mixture 
only  was  in  a  nearby  well.  Images  were  sequentially  acquired  over 
approximately  10  half-lives  (19  h,  43  min),  and  each  data  point 
was  formed  from  5  min  of  photon  accumulation. 

The  spatial  resolution  of  the  system  was  characterized  using 
a  standard  PET/SPECT  phantom  (Micro  Hot-Spot  Phantom; 
Data  Spectrum  Corp.)  filled  with  15.2  MBq  (410  pCi)  of  *®F- 
FDG;  this  phantom  is  also  commonly  used  to  evaluate  the  spatial 
resolution  of  PET  systems.  Both  the  2.4-  and  the  1.6-mm  cylin- 
dric  holes  were  used  to  form  line  profiles,  taken  from  the  fiber- 
based  images.  An  ambient-light  image  and  a  functional  Ceren¬ 
kov  luminescence  image  were  acquired.  The  exposure  time  for 
all  ambient  images  was  less  than  1  s,  and  the  functional  image 
acquisition  time  was  5  min.  The  subject  was  5  cm  from  the 
optical  system. 

Fiber-Based  CLE  System  for  Surgery  Monitoring 

Five  mice  were  injected  with  37  MBq  (1  mCi)  of  *®F-FDG  via 
the  tail  vein.  The  tracer  was  allowed  to  accumulate  for  60-70  min, 
and  the  mice  were  imaged  in  the  IVIS  system  to  verify  tumor 
uptake  of  '*F-FDG.  Tumors  were  then  imaged  with  the  IVIS 
and  fiber-based  CLE  systems  after  removal  of  the  skin  covering 
the  subcutaneous  tumor  and  after  excision  of  the  tumor.  To  sim¬ 
ulate  an  environment  mimicking  surgical  resection,  and  to  validate 
the  origin  of  the  signal,  the  excised  tumor  was  also  imaged  after 
being  placed  adjacent  to  the  surgical  cavity.  Two  images  were 
acquired  for  each  of  these  steps;  1  ambient  image  and  1  functional 
Cerenkov  luminescence  image.  The  exposure  time  for  all  ambient 
images  was  less  than  1  s,  and  the  functional  image  acquisition 
time  was  5  min. 

Statistical  Methods 

Quantitative  data  were  expressed  as  mean  ±  SD.  Means  were 
compared  using  the  Student  t  test.  A  95%  confidence  level  was 
chosen  to  determine  the  significance  between  groups,  with  P  val¬ 
ues  of  less  than  0.05  indicating  significant  differences. 


FIGURE  2.  System  sensitivity  via  sequentiai  imaging  of  3.7  MBq 
(too  pCi)  of  ■'®F-FDG  for  approximateiy  20  h.  Opticai  signais  from 
■'®F-FDG  and  controi  sampie  were  plotted. 

RESULTS 

System  Characterization 

Figure  2  shows  the  reduction  in  signal  during  '*F-FDG  [Fig.  2] 
decay  over  time  for  2  wells:  1  filled  with  ^^F-FDG  and  1 
filled  with  a  water  and  glycerol  mixture  only.  With  a  5-min 
integration  time,  a  minimum  of  approximately  45  kBq 
(1.21  pCi)  of  activity  can  he  identified  as  different  from 
the  control  well  containing  the  water  and  glycerol  solution 
(signal-to-noise  ratio  >  1).  There  was  a  decrease  in  signal 
in  the  control  well  due  to  stray  y-photons  from  the  ^*F-FDG 
well  interacting  with  the  optical  fiber  and  creating  scintil¬ 
lation  background  light  in  the  image. 

Visually,  the  photographic  and  Cerenkov  images  taken 
using  the  standard  IVIS  imaging  system  (Supplemental  Fig. 

2)  are  similar  to  the  photographic  and  Cerenkov  images 
taken  with  the  fiber-based  system  (Fig.  3A).  The  line  pro-  [Fig.  3[ 
files  demonstrate  a  high  correlation  between  the  Cerenkov 
and  photographic  line  profiles.  For  the  2.4-mm  holes  (Fig. 

3B),  the  line  profile  limits  depicted  by  the  large  green  stars 
in  Figure  3A  had  a  Pearson  correlation  coefficient  of  0.78 
(significance,  P  <  le— 5).  For  the  1.6-mm  holes  (Supple¬ 
mental  Fig.  3),  the  line  profile  limits  depicted  by  the  smaller 
red  stars  had  a  Pearson  correlation  coefficient  of  0.71  (sig¬ 
nificance,  P  <  le— 5).  The  peaks  were  also  in  near-identical 
locations;  the  difference  in  peak  locations  between  the  holes 
for  the  line  profiles  of  the  photographic  and  Cerenkov-emis- 
sion  images  was  3.8%,  a  difference  of  only  a  few  pixels. 


System  Demonstration 

The  images  comparing  both  systems  before  removal  of 
the  tumor  for  mouse  1  are  shown  in  Figure  4  (IVIS  system  [Fig.  4] 
[Fig.  4A]  and  fiber-based  CLE  system  [Fig.  4B]).  Figure  5A  [Fig.  5[ 
shows  the  images  produced  by  the  IVIS  system,  and  Figure 
5B  shows  images  from  the  fiber-based  CLE  system  after 
removal  of  the  tumor  in  mouse  1 . 

Comparing  images  in  Figure  5B  quantitatively,  we  de¬ 
termined  the  tumor-to-background  ratio  for  the  Cerenkov 
signal  by  computing  the  ratio  of  the  median  value  in  the 
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FIGURE  3.  Characterization  of  fiberscopic 
system  spatiai  resoiution.  (A)  Respective 
ambient  (ieft)  and  Cerenkov  images  (right) 
of  PET/SPECT  phantom  (4.3-cm  inner  diam¬ 
eter).  (B)  Quantitative  iine  profiies  of  ambient 
and  Cerenkov  images  for  2.4-mm  hoies;  iine 
sampied  is  indicated  by  green  stars  in  A. 
a.u.  =  arbitrary  units. 


region  of  interest  encircling  the  excised  tumor  to  the  me¬ 
dian  value  in  the  cleared  tumor  cavity.  This  ratio  was  1.28 
(for  comparison,  the  tumor-to-background  ratio  from  the 
IVIS  system  was  1.16  for  the  3-min  scan).  For  mice  2,  3, 
4,  and  5,  the  tumor-to-background  ratios  for  the  removed 
tumor  were  1.41,  1.21,  1.02,  and  1.17,  respectively.  Tumor 
tissue  light  emission  was  significantly  higher  than  the  exposed 
cavity  for  all  mice  (Student  t  test  for  paired  samples,  P  <  0.05 
for  all).  Much  of  the  low  tumor-to-background  value  in  mouse 
4  can  be  explained  by  the  reflection  of  the  light  emitted  from 
the  tumor  by  the  tumor  cavity,  which  was  directly  adjacent. 
Residual  tumor  tissue  after  surgery  is  also  possible. 

DISCUSSION 

In  this  study,  we  demonstrated  the  methodology  and 
performance  of  an  optical  fiber  system  built  to  image  tumor 
margins  during  surgical  resection.  Cerenkov  imaging  has 
been  proposed  for  small-animal  drug  discovery  and  surgical 
guidance,  yet  previous  studies  have  used  an  imaging  box 
with  a  lens-mounted  imaging  system;  this  setup  is  not  realistic 
for  subjects  that  are  larger  (e.g.,  large  animals  and  humans). 
To  our  Icnowlcdgc,  this  was  the  first  study  that  investigated 

a  Cerenkov  endoscope  for  surgical  resection.  We  have  shown 
that  this  system  is  capable  of  imaging  1.2-mm  structures, 
using  phantom  experiments.  We  determined  the  sensitivity  of 
the  system  for  '^F-FDG  (—45  kBq  [1.21  p,Ci]/300  p.L).  We 
then  demonstrated  the  ability  to  visualize  the  accumulation  of 
'*^F-FDG  in  a  tumor  using  this  system  and  could  track  the 
contrast  as  the  tumor  was  excised  from  the  mouse. 

A  benefit  of  a  Cerenkov  endoscope  is  that  it  provides 
higher  image  spatial  resolution  than  PET  and  SPECT 
cameras,  because  the  optical  photons  emitted  from  the 
tissue  surface  are  detected  with  optical  lenses  and  a  high- 
resolution  CCD  chip.  Cho  et  al.  demonstrated  a  limit  of 
350  p-m  (in  full  width  at  half  maximum)  with  a  system 
designed  for  visualizing  a  microfluidic  chip  setup  (-3^). 
This  resolution  is  affected  by  the  imaging  system,  the 
range  of  the  (3-particles  in  the  tissue,  and  tissue  optical 
scatter.  Levin  and  Floffman  showed  that  the  mean  (3 -particle 
track  of  ^®E-EDG  had  a  full  width  at  half  maximum  value 
of  102  p,m  in  water  (-34),  which  determines  the  physical 
resolution  limit  of  a  Cerenkov  scope  when  imaging  '*E- 


FDG.  Improved  resolution  may  be  achieved  using  a  radio¬ 
isotope  with  a  (3-particle  decay  of  lower  energy,  such  as 
^^'I,  although  this  would  result  in  lower  sensitivity  be¬ 
cause  there  is  an  inverse  relationship  between  resolution 
and  sensitivity.  We  do  note  some  inconsistency  in  the 
image  of  the  1.2-mm  holes,  but  many  of  these  holes  are 
clearly  visualized.  The  high  spatial  resolution  of  the  sys¬ 
tem  requires  accurate  focus  of  the  optics  both  at  the  distal 
fiber-lens  junction  and  at  the  proximal  fiber-camera  junc¬ 
tion — an  optical  focus  that  will  be  improved  in  our  further 
studies.  In  the  case  of  our  phantom,  the  resolution  of  the 
Cerenkov  luminescence  images  was  also  degraded  by 
reflections  within  the  translucent  phantom  and  the  sig- 
nal-to-noise  ratio  of  the  single-photon-counting  camera. 


FIGURE  4.  Mouse  1  bearing  C6  glioma  after  tail-vein  administra¬ 
tion  of  37  MBq  (1  mCi)  of  ''®F-FDG.  (A)  Mouse  was  imaged  by 
commercially  available  optical  IVIS  system,  and  images  were  com¬ 
pared  with  those  from  prototype  fiber-based  system  (B).  Tumor 
tissues  are  outlined  by  red  lines.  Ambient-light  images  are  on  left, 
luminescent  images  are  in  middle,  and  fused  images  are  on  right. 
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FIGURE  5.  Mouse  1  was  imaged  by  IVIS  optical  system  (A)  and 
fiber-based  system  (B)  after  surgery  to  remove  tumor  tissues. 
Ambient-light  images  are  on  left,  luminescent  images  are  in  middle, 
and  fused  images  are  on  right. 


endoscopic  or  laparoscopic  applications.  However,  smaller 
custom  optics  could  be  designed  to  enable  imaging  through 
smaller  apertures.  Meanwhile,  optimized  optic  lenses  and 
fibers  can  provide  better  sensitivity  and  shorter  integration 
time  in  upcoming  systems.  We  leave  this  investigation  for 
future  studies  using  more  customized  optical  systems. 

The  main  advantage  of  this  system  over  other  optical 
techniques,  such  as  fluorescence  imaging,  is  that  it  is  able 
to  use  standard  PET  tracers,  such  as  ^^F-FDG,  which  are 
already  available  in  the  clinic.  Many  novel  tracers  such  as 
3'-deoxy-3'-^*F-fluorothymidine,  engineered  proteins,  and 
labeled  arginine-glycine-aspartic  acid  peptides  are  in  the 
pipeline  of  clinical  trials  required  by  the  Food  and  Drug 
Administration  because  of  their  utility  in  PET.  This  is  a  sig¬ 
nificant  advantage  for  this  technique,  because  much  effort 
and  innovation  is  being  dedicated  to  PET  tracers.  In  con¬ 
trast,  because  of  the  lack  of  a  standardized  imaging  device 
for  optical  imaging,  optical  contrast  agents  are  not  being 
sufficiently  developed  to  keep  pace  with  radiotracers.  Thus, 
the  main  advantage  of  a  Cerenkov  emission  surgical  scope 
is  that  it  is  readily  clinically  translatable.  The  combination 
of  PET  and  Cerenkov  imaging  provides  an  opportunity  to 
use  an  identical  signal  origin,  and  perhaps  an  identical  in¬ 
jection,  to  localize  tumors  for  excision. 


which  was  not  entirely  optimized.  The  signal-to-noise  ra¬ 
tio  of  our  ambient-light  image  was  also  affected  by 
a  slightly  imperfect  focus  and  the  suboptimal  settings  of 
the  illumination  source  and  camera.  Thus,  although  we 
demonstrated  the  ability  to  visualize  many  of  these  1.2- 
mm  cylindric  holes  in  a  phantom  with  this  endoscope- 
based  system,  we  expect  to  be  able  to  image  submillimeter 
lesions  with  a  more  optimized  system. 

The  sensitivity  limit  of  this  system  was  45  kBq  (1.2 
p.Ci)/300  p-L  with  '*F-FDG.  To  put  this  in  a  practical  per¬ 
spective,  if  this  system  were  used  for  head  and  neck  tumors 
of  the  oropharynx  (with  a  standardized  uptake  value  of 
5.17),  2.1  GBq  (54  mCi)  would  have  to  be  injected  intra¬ 
venously  into  a  70-kg  patient  to  have  enough  activity  to 
visualize  a  300-p.g  tumor  (35)  in  a  5-min  scan;  a  more 
sensitive  scope  would  enable  the  identification  of  smaller 
structures  or  similar  structures  with  a  reduction  in  ionizing 
dose  to  the  patient  and  clinical  staff.  Increased  sensitivity 
may  be  realized  with  more  sensitive  optics  (such  as  an 
F/0.95  lens)  or  fiber  optic  glass,  which  transmits  farther 
into  the  violet  and  ultraviolet  emissions.  The  fiberglass  used 
in  this  experiment  was  Schott-75  glass,  which  transmits 
40%  light  at  500  nm;  thus,  it  is  not  optimized  for  detecting 
the  short-wavelength-dominant  Cerenkov  emission.  In  ad¬ 
dition,  y-photons  emitted  from  the  radiotracer  scintillate 
impurities  in  the  glass,  resulting  in  undesired  background 
noise.  The  use  of  fused  silica,  which  is  much  more  sensitive 
to  the  ultraviolet  and  violet  emissions  and  has  fewer  impu¬ 
rities,  would  improve  sensitivity. 

As  a  prototype  device,  the  off-the-shelf  optical  lens 
on  this  system  was  too  large  (3-cm  diameter)  for  many 


CONCLUSION 

This  proof-of-concept  study  demonstrated  the  feasibility 
of  using  fiber-based  CLE  for  the  detection  of  tumor  tissue 
in  vivo  and  demonstrated  its  potential  use  for  image-guided 
surgery.  With  further  improvement  in  imaging  sensitivity  and 
resolution  of  the  current  system,  it  is  expected  that  Cerenkov 
imaging  might  soon  be  translated  into  clinical  applications. 
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